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We present wide- field near-infrared images of the DR21/W75 high-mass star form- 
ing region, obtained with the Wide Field Camera, WFCAM, on the United Kingdom 
Infrared Telescope. Broad-band JHK and narrow-band H2 l-OS(l) images are com- 
pared to archival mid-IR images from the Spitzer Space Telescope, and 850iJ.m dust- 
continuum maps obtained with the Submillimeter Common User Bolometer Array 
Q . (SCUBA). Together these data give a complete picture of dynamic star formation 

' across this extensive region, which includes at least four separate star forming sites 

^ \ in various stages of evolution. The H2 data reveal knots and bow shocks associated 

. with more than 50 individual flows. Most are well collimated, and at least five qualify 

as parsec-scale flows. Most appear to be driven by embedded, low-mass protostars. 
The orientations of the outflows, particularly from the few higher-mass sources in 
the region (DR21, DR21(0H), W75N and ERO 1), show some degree of order, be- 
5— i ' ing preferentially orientated roughly orthogonal to the chain of dusty cores that runs 

north-south through DR21. Clustering may inhibit disk accretion and therefore the 
production of outflows; we certainly do not see enhanced outflow activity from clusters 
of protostars. Finally, although the low-mass protostellar outflows are abundant and 
widely distributed, the current generation does not provide sufficient momentum and 
kinetic energy to account for the observed turbulent motions in the DR21/W75 giant 
molecular clouds. Rather, multiple epochs of outflow activity are required over the 
million-year timescale for turbulent decay. 
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outflows - ISM: kinematics and dynamics - ISM: DR21 - ISM: DR21(0H) - ISM: 
W75 
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1 INTRODUCTION 

Star formation studies have in recent years benefited from 
a fresh, "large scale" view, via observations of Giant Molec- 
ular Clouds (GMCs) and large clusters of clumps and cores 
at infrared and (sub)millimetre wavelengths, and through 
studies of the interplay between the thousands of evolving 
young sources within each region. 

Infrared survey instrumentation that ofi'ers an unpar- 
alleled field of view at high spatial resolution is now, or 
will soon be, available at modest-sized telescopes across the 



globe. One such instrument is the near-infrared (near-IR) 
Wide Field Camera, WFCAM, at the U.K. Inft-ared Tele- 
scope (UKIRT). WFCAM provides deep imaging of degree- 
sized regions in broad and narrow-band filters in just a few 
hours. When combined with photometric data at longer 
wavelengths, e.g. from the Submillimeter Common User 
Bolometer Array (SCUBA - soon to be replaced by SCUBA- 
2) and the Spitzer Space Telescope, obtaining a complete 
picture of evolving star formation (stellar populations, the 
IMF, the clump mass function, etc.) becomes possible. Com- 
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plementary wide-field optical images of Herbig-Haro (HH) 
jets, and future large-scale (sub)millimetre surveys of molec- 
ular "CO" outflows, may also be used to examine the dy- 
namic processes that shape - through the injection of tur- 
bulence or triggered cloud collapse - star formation across 
each CMC. 

One region that has already benefited from wide-field 
observations with Spitzer, SCUBA and now WFCAM is 
the DR21/W75S, DR21(OH) and W75 N complex of mas- 
sive star forming cores in Cygnus X jWesterhoud llQSSt 
iDownes fc Rinehartlll96^ . for a more recent overview of the 
Cygnus X complex of GMCs see Schneider et al. 2006). 
The D R21 core con t ains two well-studied cometary Hll 
regio ns ijHarri j HqtS IPickel et al.lll98fit ICveanowski et alJ 
l200,?fl and an H2O maser ijGen^elfcPowne 
while DR21(OH), a dense core roughly 3' to the north 
of DR2 1, is abundant with me thanol, water and OH 
masers iLiechti fc WalmslevUlQQ^ and references therein). 
DR21(0H) is undetected in radio continuum emission and 
so is probably younger than DR21. W75N, roughly 20' to 
the north of DR21, comprises a group of embedded infrared 
sources (M_o_orc_ot_al. 1991b; Porsi, Tapia & S mith 20o3), 
Hll regions and compact radio sources (e.g. iHaschick et alJ 
Il98ll: iHunter et al. | ll994r). as weU as H9 O and OH masers 
jHunter et al.lll994l:nbrrelles et alJll99'it . 

E arly work by iDickel et alj ([l97^ and 'Fisc her et alJ 
^1985^ showed that DR21 and W75N coincide with two, 
possibly interacting molecular cores along a north-south 
molecular ridge extending over at least half a degree. These 
cores, evide nt in the more recen t 85 lira dust contin- 
uum maps of IVallee fc Fieed J2006D and lGibb et"ai] J2007D . 
have been resolved into smaller clumps by various groups 
iWilson fc Mauersbergerll 990l:lMangum. Wo oten fc Mundvl 
Il992l: IChandler. Gear fc Chinil Il993al : IShepherdI I2OOID . 
Spitzer images at 5.8 and 8.0 fim reveal streamers of 
emitting material that seem to radiate away from the 
brighter infrared sou rces located along the molecular ridge 
iMarston et alJl2'004^ . features possibly produced by ioniza- 
tion and ablation of the less dense material on either side 
of the ridge. .Marston e t al. ( 2004) also identify a number of 
"extremely red objects" , or EROs, along the ridge. These are 
luminous, highly-reddened objects that may be very young, 
massive, early B-type sources. 

In terms of outflows, DR21 and W75N have been stud- 
ied independently by a number of groups. DR21 drives 
one of the most spectacular flow s known iCarden et alJ 
Il991 b': 'Gard en fc Carlstromlll992ll . T he bipolar outflow is 
extre mely b right in H2 line emission JCarden et alJll991al: 
iDavis fc SmithI Il996l : ISmith. Eisloffel fc Davis! Il998^ and 
produces methan ol abundance enhancement s and masers 
through shocks iLiechti fc Walmsle^ [l993). The domi- 
nant outflow from W75N is somewhat fainter in H2, al- 
though it appears to be a good example of a wind-driven 
flow, in which sweeping bow shocks e ntrain ambient gas 
to form a massive molecular outflow jDavis et alJ Il998al : 
iDavis. Smith fc Moriartv-Schievenlll998b^ . Earlv CO maps 
showing the orientation and bipolarity of the W75N out- 
flow were presented by Fischer et al. (198^; more de- 
t ail ed maps i n CO and millimetre continuum em ission 
bv IShepherdI j200lD. Ishepherd. Testi fc Stark! (|200i) and 
IShepher^^^r^fcTtes^ i2004) illustrate the complexitv of 
the region. 



In this paper we discuss WFCAM data obtained in H2 
l-OS(l) emission and in the J, H and K bands. The images 
cover a 0.8° x 0.8° "tile" (see Section 2) centred close to 
DR21. The broad-b and photo metry are discussed in detail in 
a companion paper iKumar e t al. 2006, hereafter Paper II); 
in this work we focus on the H2 emission-line features associ- 
ated with the many outflows along the DR21/W75N ridge. A 
comparis on is a l so ma de wi th SCUBA 850 lira data (see also 
.Vallee fc Fieed ll200e!) and ICibb et al.l (|20o3) for a similar 
map), and archival Spitzer IRAC images at 3.6 ura, 4.5 ura, 
5;8um and 8.0 ^ra (originally presented by iMarston et alJ 
I2TOJ. The combined SPITZER and SCUBA maps suggest 
that stars are forming within a flattened north-south flla- 
ment viewed almost edge on. One of the goals of the WF- 
CAM observations was to search for outflows ejected east- 
west, i.e. perpendicular to this dense filament. We discuss 
the effect these outflows might have on the general CMC, 
and look for examples of collimated jets from the more mas- 
sive YSOs. 

There has been some disagreement on the distance 
to DR21 and W75N, although both arc thought t o be in 
the range 1.5-3.0 kpc jC enzel fc Downcs 1977; F ische^ t alJ 
Il985ll0denwald fc Schwartzll993l : ISchneider et al.l2006h . In 
Paper II we flnd that the upper limit gives a more realistic 
population of massive stars, so we adopt 3 kpc as a general 
distance in this paper. 



2 OBSERVATIONS 

2.1 WFCAM and Spitzer imaging 

Wide- field images of the DR21/W75 region were ob- 
tained during service observing at the United Kingdom 
Infrared Telescope (UKIRT) on 6 May, 11 May and 9 
June 2005, using t he near-IR wide-field camera WFCAM 
iCasali et al.l I2OO6I) . The camera employs four Rockwell 
Hawaii-II (HgCdTe 2048x2048) arrays spaced by 94% in the 
focal plane. The pixel scale measures 0.40". 

To fill in the gaps between the arrays, and thereby ob- 
serve a contiguous square field on the sky covering 0.75 
square degrees - what is referred to as a WFCAM "tile" 
- observations at four positions are required. At each posi- 
tion, to correct for bad pixels and array artifacts, a five-point 
jitter pattern was executed (with offsets of 3.2" or 6.4"); to 
fully sample the seeing, at each jitter position, a 2x2 micro- 
stepped pattern was also used, with the array shifted by half 
a pixel. In this way 20 frames were obtained at each of the 
four positions in the tile. 

Data through Mauna Kea Consortium broad-band J, H 
and K filters, and a narrow-band H2 l-OS( l) filter (A = 
2.121 /xm, &\ = 0.021 /im), were obtained jHewett et alJ 
200§). Exposure times of 5 sec and 40 sec were used with 
the broad and narrow-band filters, respectively. A correlated 
double sampling (CDS) readout mode was used with the 
broad-band filters, while a non-destructive readout (NDR) 
mode was employed with the longer-exposure H2 images. 
The latter yields a lower read noise (~20e- as compared to 
~30e- with CDS) at the expense of greater amplifler glow 
and a higher risk of electronic pickup. With the 5-point jitter 
pattern and 2x2 micro-stepping in each broad-band filter, 
the total on-source/per-pixel integration time was 100 sec 
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per filter. In H2 the whole 20-frame sequence was repeated; 
hence, the total per-pixel exposure time was 1600 sec. 

The data were reduced by the Cambridge Astronomical 
Survey Unit (CASU), which is a fore-runner to the VISTA 
Data Flow System Project (VDFS), a collaborative effort be- 
tween the Universities of London, Cambridge and Edinburgh 
in the U.K. CASU is responsible for the design and imple- 
mentation of the data processing pipeline used prior to the 
archiving and subsequent distribution of all WF CAM data. 
The r educ tion steps are desc ribed in detail by iDve et all 
J200d) and llrwin et alJ i2006i) . Briefly, dark frames secured 
at the beginning of the night were subtracted from each 
target frame to remove bad pixels, amplifier glow, and re- 
set anomalies from the raw frames. Twilight flats were then 
used to flat-field the data. The micro-stepped frames were 
interleaved and the offset frames combined by weighted av- 
eraging (using a confidence map derived largely from the 
flat-field frames) to produce "leavestack" frames. 

Astrometric calibration was achieved through a cubic 
radial fit and a six-parameter linear transformat ion to ob- 
jccts found in the 2MASS point source catalogue jPve et alJ 
b006). The transformation takes into account translation, 
scaling, rotation and shear. Photometric cal ibration was also 
achieved through the use of 2MASS data jPve et alJl2'ooi 
iHewett e t al. 2006). 

AU reduced WFCAM data are archived at the WF- 
CAM Science Archive (WSA), which is part of the Wide 
Field Astronomy Unit (WFAU) hosted by the Insti- 
tute for Astronomy at the Royal Observatory, Edinburgh, 
U.K. (http://surveys.roe.ac.uk/wsa/index.htmll. WSA data 
are distributed as RICE-compressed Multi-Extension Fits 
(MEF) files. The data for one flher covering the 0.8° xO. 8° 
tile appear as four leavestack files in the database, one MEF 
file for each of the four offsets on sky needed to cover the 
tile. Each leavestack file therefore contains four individual 
images, one for each of the four WFCAM cameras. 

Having retrieved the reduced data from the WFCAM 
Science Archive (WSA), Starlink software were used to re- 
move residual large-scale gradients across each of the 16 
images in each tile. The Starlink:KAPPA command SUR- 
FIT was used to flt a coarse surface to the background in 
each frame; the four surface fit images in each corner of 
the full tile were median- averaged and used to flat-fleld the 
four associated target frames. The frames were subsequently 
mosaicked together using the KAPPA:WCSALIGN routine 
applied to the world coordinates assigned to each frame by 
the CASU pipeline, together with the Starlink:CCDPACK 
mosaicking routine MAKEMOS. These gave one complete 
tile image in each filter. 

The Spitzer Space Telescope images discussed in this 
paper were retrieved from the Spitzer archive. Our process- 
ing of these data is described in Paper II. The observations 
were obt ained with the IRAC camera, which is described by 
^zi^^^al . 12004' ). The dat a also appear in recen t pape rs by 
Marston 618111200 4) and Persi. Tapia fc SmithI (120061) . al- 
though |M^sto^^L.aLt (|200j) only give preliminary results, 
while iPersi. Tapia fc SmithI ll2006f) focus only on a small 
-2'x2' field centred on W75N. 



2.2 SCUBA observations 

The sub-millimetre data were obtained wit h the SCUBA 
bolometer array camera (iHolland et alJll999D on the James 
Clerk Maxwell Telescope (JCMT) on Manna Kea, Hawaii 
on four partial nights from May 26 to May 29, 2001 in 
marginal sub-millimetre conditions; the precipitable water 
vapor (PWV) was ~1.3-2 mm (optical depth at 850 /^m 
^0.25-0.36), although it was fairly stable each night. Due 
to the high PWV content in the atmosphere only the 
850 ^tm data were usable. The standard "scan-mapping" 
mode jjenness. Lightfoot fc Holland! Il99d) was employed 
with six different chop throws: 30", 44", and 68" in Right As- 
cension and Declination respectively, while scanning in the 
Nasmyth coordinate frame. Each chop throw was repeated 
2-4 times depending on airmass and sky opacity, the goal 
being to reach the same noise level in each dual beam map. 
The rms noise level in the final co-added dual beam maps 
was ~0. 09-0.1 Jy beam~^. The pointing was checked fre- 
quently on the point source MWC 349, and each dual beam 
map was corrected for any pointing drift between pointing 
observations. Based on the pointing data the positional ac- 
curacy in the final 850 /im image is estimated to be ^1". 
The maps were calibrated with scan maps of Uranus as the 
primary flux calibrator and CRL 2688 as the secondary cali- 
brator. The flux calibration is estimated to be good to within 
5% with respect to Uranus. The telescope half power beam 
width was measured to be 15.2" from scan maps of Uranus. 

The maps were reduced in a standard way us- 
ing the Starlink reduction pac kages SURF and Kappa 
iSandell. Jessop fc Jennes^l200lll . except that the individ- 
ual dual beam maps where created by weighting the noise 
in each individual bolometer using the SURF task setbolwt, 
which reduced the noise level in the maps by ~5-10% . We 
then estimated the rms noise level using emission free areas 
in the maps and co-added the maps by accounting for the dif- 
ference in noise between individual maps for the same chop 
throw and chop position angle. The flnal calibrated 850 ^im 
map has a noise level of 50-60 Jy beam~^ over the whole 
area. For further analysis we divided the map into three 
areas (south, middle, and north), removed the error beam 
using the MIRIAD task CLEAN dSauh. Teuben fc Wright 
1995j and restored the maps to 14.0" resolution. Below we 
present only photometry of the major 850 /im cores, plus 
contour maps for comparison with the near- and mid-IR 
data. A more thorough analysis of the 850 fj,m observations 
is beyond the scope of this paper, and will be presented in 
a later work. 



3 RESULTS 

3.1 WFCAM imaging of DR21/W75 

Our H2 l-OS(l) (4- continuum) tile of the DR21/W75 region 
is presented in Fig.l (which shows the full field observed). 
Regions where H2 line-emission features were identified have 
been outlined. Although somewhat clustered around DR21 
and W75N, these generally follow the ridge of obscuration 
that links DR21 in the centre with W75N in the north- 
northwest. There is also a cluster of features to the west- 
south-west of DR21, which may be unrelated to the main 
high-mass star forming ridge. 
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Figure 1. A near-IR narrow-band H2 2.122 image of tlie DR21 and W75N regions. The figure shows a full WFCAM tile, covering 
an area of roughly 0.8x0.8 degrees. Areas were H2 line emission is identified are outlined; the regions A, B, C, D and E are discussed 
separately in Sect. 4. 



Colour images of the regions around DR21 and W75N 
are shown in Figs. |21 and |3 J, K and H2 mosaics were 
scaled by exposure time, filter bandpass and filter trans- 
mission before being registered and combined. In the re- 
sulting figures embedded sources can be easily distinguished 
from pure line-emission features, i.e. the H2 knots and bow 
shocks in jets and outflows. The former appear yellow while 
the latter are bright red in colour. Some of the yellow fea- 
tures in the DR21 lobes are due to scattered continuum, 
while some may result from "boosted" K-band emission due 
perhaps to enhanced vibrational excitation by fluorescence 



iFernandez. Brand fc Burton|[l997l) . Many of the features 
marked in Fig. 1 are seen in more detail here; indeed, these 
colour images were used to search for H2 features across 
this extensive region. In Fig.|5|the molecular ridge that runs 
through the middle of DR21 is evident as a dark lane, largely 
devoid of point sources. We mark on Fig. |21the "Extremel y 
Red Objects", or EROs, identifled bv lMarston et alJ J2004) . 
as well as the location of the DR21 A-B-C and D Hll regions 
JCveanowski et al.ll2003D . and the DR21(OH) cluster of wa- 
ter, OH and methanol masers (iLiechti fc Walmslev'l 997i). In 
Fig.Olthe Hii regions A, B and C in W75N (Haschick et alJ 
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Figure 2. A colour image of tlie DR21 region composed of observations in J (blue), K (green) and H2 (red). With this combination 
H2 line emission regions appear red, while embedded or background stars appear yellow. The star and triangle mark the DR21 A-B- 
C and DR21 D cometary Hll regions; the cro ss indicates the locat ion of the DR21(OH) cluster of water, methanol and OH masers 
jLiechti k Walmslevll99^ . EROs identified bv lMarston et akl i2004) are labeled; the infrared "DR,21-IRS" sources discussed in the text 
are numbered. 



Il98ll) are labeled. In both figures K-band sources with bright 
8.0 /xm counterparts are numbered; these are probably the 
most massive, embedded sources in the region. 

DR21 is famous for its Hll regions. It comprises four 
ma in radio emission peaks, originally labeled A, B, C and 
D llHarrislll97i iRoelfsema. Goss fc Geballel|l989^■ None of 
the r adio peaks are detected in the near-IR 1 Davis fc SmithI 
Il99(fl , although the nort hern, compac t sourc e D is observed 
at mid-IR wavelengths. ISmith et alJ ll2005l) point to this 
source as being the dominant source of infrared luminos- 
ity in DR21. They find that its total luminosity exceeds 
that of an 08 ZAMS star and attribute the excess to ac- 
cretion. They estimate an extinction of ~ 175 mag to- 
ward source D, and suggest that because it is still accret- 
ing, it could (at least in part) be responsible for the large- 
scale DR21 outflow. However, at radio wavelengths Hll re- 
gion D is c ometary in appearance, opening outward toward 
the north iCveanowski et al]|2003l) . One could argue that 
a north-south flow is more likely. Either way, there are no 



H2 knots or jet-like features close to this source which one 
might confidently associate with the Hll region. 

Situated ~20" south of source D, radio sources A, B 
and C appear to be peaks in a larg er "bow shock" or "cham - 
pagne flow" style cometary nebula iCveanowski et al.l2003l) . 
This nebula lies ^-^15" to the east of the cluster of near-lR 
sources in DR21, within the dark lane that separates the 
two lobes of the DR21 outflow. It is coincident with a dense 
molecular core (see Section 3.3), although it has no mid-IR 
counterpart. Overall, the cometary A-B-C nebula opens in 
the direction of the eastern lobe of the main DR21 outflow. 
This might suggest a causal relationship between the flow 
and the Hll region. However, the more coUimated, western 
lobe of the DR21 outflow is difficult to reconcile with the 
sharp, western edge of the cometary Hll region. The power- 
ful flow would surely disrupt this Hll region/molecular cloud 
boundary, if the outflow central engine was situated within 
the cometary nebula. 

The complex of masers and mm-wave contin- 
uum peaks known as DR21(OH) are found about 3' 
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Figure 3. Same as Fig.|21 except for the region around W75N. Tlie brigiit infrared sources referred to as W75N-IRS 1, etc. in tlie text 
are numbered. Tlie letters A, B and C refer to the main Hn regions. 



to the north of DR21 (A-B-C-D). DR2 1(OH) coin- 
cides with a chain of four water ma sers iForster et alJ 
Il978l: iManeum. Wooten fc Mundv^ll992^ spread northeast- 
southwest over about 20". Orthogonal to the water masers, 
about a dozen methanol masers dchneate the two lobes 
of the DR21(OH) outflow (iPlambcck fc Menteiil Il990l: 
iKoean fc Slvshl 119981) . The water masers likely trace em- 
bedded young stars and ultra-compact Hll regions, includ- 
ing in this case the outflow central engine, while the class I 
methanol masers probably trace the interaction between the 
outflow and the ambient medium. H2 line emission is also 
excited in such outflow-ambient impact zones, although we 
fail to detect clear evidence for H2 emission in our WFCAM 
data. The absence of H2 emission in the DR21(0H) outflow 
is likely due to extinction. 

W75N, shown in Fig. |21 comprises three regions 
of radio emission, W 75N(A), W75N(B) and W75N(C) 
I'Haschic k et al.|[l98 j) . A and B are coincident with patches 
of bright near-IR nebulosity; the more evolved source, 
W75N(A), is associated with a 2 /xm point source. The over- 
all cluster harbours at least 30 IR sources with infrared 
excesses, probably due to circums tellar disks (iMoore et alJ 
ll991bl : IPersi. Tapia fc Smithll200d) . 



Collectively, W75N is associated w i th a ~ 10^ M© 
dense molecular core llMoore et al.l liggiaTl . which 
IShepherd. Testi fc StarS ||200^ resolve into nine dis- 
crete millimetre continuum peaks spread over an area 
of about an arcminute. The brightest and most massive, 
MM 1, peaks a few arcseconds to the northwest of IRS 1 
and is associated with Hll region W75N(B); W75N(A) is 
similarly associated with a continuum peak, MM 5. MM 1 
contains at least four Hll regions (a fifth is associated with 
MM 5). Thre e of these were de tected at 1.3 mm with the 
VLA by Torrefies et alj Jl997l) . Of these, VLA 1 drives 
a radio jet in the direction of the large-scale molecular 
outflow, althoug h all three VLA sources may drive inde- 
pendent flows dShepherdl I2OOII: ISheoherd. Testi fc Stark! 
i2003). Al l three sources are a ssociated with water maser 
emission llTorrelles et alJll997^ : VLA 1 is also surrounded 
by a compact cluster of OH rn asers iBaart et al.l Il986l: 
iHutawarakorn. Cohen fc BrebneH l2002h . Each source is 
presumably in the protostellar phase, their radio fluxes 
being consistent with ea r ly B- type ZAMS stars. Indeed, 
IShepherd. Kurtz fc Testil 11200411 assign spectral types of 
B0-B2 to ah flve Hii regions, W75N(A) being perhaps the 
most evolved at BO. 5. Our H2 observations of W75N reveal 
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Figure 4. Left: WFCAM H2 2.12/im image of DR21 with Spitzer IRAC 8.0/im contours ovcrlayed. The contours pick out only the 
brightest 8.0/im peaks; lower contours associated with diffuse PAH emission are not plotted. Contour levels are 3.4, 5.1, 6.8, 10.2 mjy 
arcsec"-^ (black) and 17.0, 30.6, 85.0, 165.0 mJy arcsec"^ (white); inset, the same contours are plotted but with thin and thick lines. 
DR21-IRS sources are numbered. Right: Spitzer IRAC 8.0 ^m image with SCUBA 850 ^m contours overlayed. Contour levels are 0.125, 
0.25, 0.5, 1, 2, 4, 8 and 16 Jy/beam; the beam size measures 14". 



not only the well-known H2 features associated with W75N 
itself, but also collimated jets to the north and west. 

Lastly, roughly 15'-20' to the west of DR21 there is a 
further region of active star formation that has to date at- 
tracted little attention. The re gion is evident in the Spitzer 
data of iMarston et al.1 i2004) . although they do not dis- 
cuss these features in any detail. A search of the SIMBAD^ 
database suggests that the region is situated on the eastern 
edge of the small L906 dark cloud. We find an abundance of 



H2 hue emission knots and jets, plus regions of bright, neb- 
ulous, mid-IR emission and a number of highly reddened 
sources (Paper II). We use the suffix L906E for the sources 
in this region. 



^ The SIMBAD database is operated at CDS, Strasbourg, 
France. 



Notably, no H2 flows were found in the Diamond Ring 
region (Fig.0, which is consistent with this area being more 
evolved than DR21, W75N and L906E. 
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3.2 Comparison with Spitzer data 

3.2.1 The young stars 

In Figs. 0J[^ we compare our WFCAM observations 
of DR21/W75 with the ar chival Spitzer 8.0 ob- 
servations first discussed by iMarston et al.l ll2004f) and 
iPersi. Tapia fc SmithI l|2006^ . These longer- wavelength data 
are potentially very useful for tracing luminous, embedded 
sources. Photometry of the brighter 8.0 ^im point sources is 
given in Table Q 

The EROs identified bv lMarston et al.l J2004) in DR2f 
are bright, reddened 8 pm sources with spectral indices {a — 
dlogXFx / dlogX) between 2.2 nm and 8.0 fim in the range 1.8 
to 3.6, mu ch larger th an the a > 0.3 lower limit for Class I 
protostars jLadallQSTT) . ERO 1 lies at the western edge of an 



"infrared dark cloud" that extends eastward for at least 2' 
iMarston et al...2004l ; the cloud is evident as a region devoid 
of stars in Fig. |5| In the 850 ^m data described below the 
cloud is seen as two peaks (FIR 5 and FIR 5b) enveloped in 
an easterly elongation of the dust emission. 

The other extremely red objects, ERO 2, ERO 3 and 
ERO 4, all lie along the north-south chain of dense cores 
that runs through the centre of DR21. In Table^we list the 
Spitzer magnitudes of these and many other bright 8.0 fj.m 
sources in DR21/W75: in Fig. (around DR21), Fig.|Hl(near 
to W75N) and Fig. (around L906E), these sources are la- 
beled DR21-IRS 1, DR21-IRS 2, etc., W75N-IRS 1, W75N- 
IRS 2, etc, and L906E-IRS 1, L906E-IRS 2, etc, respec- 
tively (the same sources are labeled in the near-IR images 
in Figs. Inland ISJ. These "IRS" sources are generally bright 
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Figure 6. WFCAM H2 2.12 /^m image of tlie eastern edge of L906 (west of DR21). Spitzer IRAC 8.0/im contours are overlayed. Contour 
levels are 3.4, 5.1, 6.8, 10.2 mjy arcsec"^ (black) and 17.0, 30.6, 85.0, 165.0 mjy arcsec"^ (white). The L906E-IRS sources are numbered. 



8.0 /im point sources or luminous, reddened K-band sources 
associated with bright 8.0 emission (surface brightness 
>7 mJy arcsec"^). All appear as peaks or discrete sources 
in Figs. I4I6I Spectral indices, a, derived from linear fits to 
the Spitzer photometry between 3.6 /im and 8.0 /im (as de- 
scribed in Paper II), are also given in Table [T] We do not 
correct for extinction since this is not know a priori for each 
source, although it is worth noting that, at a distance of 
3 kpc the interstellar ex tinction to DR21/W75 may be as 
large as ~ 8.8 mag Jjoshlliooi). Extinction will have 
the effect of increasing the range of a values (Paper II). For 
targets with photometry in only two IRAC bands we do not 
derive values for a, since the 8.0 /xm magnitudes may be 
compromised by deep silicate absorption in these luminous, 
embedded sources (e.g lGibb et"ai]|2000h . 

Most of the targets in Table Q are very red. Indeed, 
when considering all of the sources, the EROs identified by 
iMarston et alJ ll2004l) do not seem to be particularly unique. 
We find many other examples of extremely red objects. 
DR21-IRS 5 and L906E-IRS 1 are particularly noteworthy, 
having a (uncorrected for extinction) of ~ 4.1. UCHll re- 
gion DR21(D) is also very red, as is DR21(OH). For the 
EROs, the spectra l indices in Tab l e U di ffer somewhat from 
those reported bv lMarston et all (|20o3), probably because 



different photometry apertures (and sky annuli) were used 
with these moderately extended sources. Our Spitzer pho- 
tometry in the W7 5N regions matches that reported by 
IPersi. Tapia fc Smith ( 2006) , who use a similar photometry 
aperture to ourselves, reasonably well. 

The mid-infrared sources in the W75N region are shown 
in Fig. 13 The region associated with W75N(A) and (B) is 
very complex. Ground-based and Spitzer mid-IR photome- 
try of the sources in the i mmediate vicinity of W75N have 
been discussed in detail bv lPersi. Tapia ilSmith (^06), al- 
though there are additional mid-IR sources to the north and 
west close to bright infrared nebulae and regions of H2 fine 
emission. Porsi, Tapia & Smith (2006) identify mid-IR coun- 
terparts to the near-IR sources W75N-IRS 1, 2 and 3, as well 
as the VLA 3 UCHll region in MM 1. They derive an extinc- 
tion of 4v ~ 90 to VLA 3. 

Much hke the IRS sources near DR21, around W75N 
and to the west, the combined WFCAM and Spitzer data re- 
veal nebulous sources which harbour a luminous O or early- 
B star within a compact cluster of less-luminous sources (e.g. 
W75N-IRS 1 to W75N-IRS 5, W75N-IRS 7 and probably 
W75N-IRS 8), as well as examples of single intermediate- 
mass sources within an almost spherical shell of reflected 
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Figure 7. Colour image constructed from WFCAM H2 data (2.12 /^m - blue) and Spitzer IRAC images (3.6 /^m - green; 4.5 /xm - red). 
EROs are labeled and DR21-IRS sources are numbered. 



near-IR light and 8.0 fim PAH emission (e.g. W75N-IRS 6, 
IRS 13 and IRS 14). 

In Fig.lHlwe show WFCAM 2.12 fim and Spitzer 8.0 fim 
data of star forming regions to the west of DR21. The Spitzer 
data reveal a chain of three spherical cores. As one moves 
westward along this chain the cores become less diffuse and 
more centrally peaked at 8.0 ^m: note that the eastern core 
in Fig. |S| has no 8.0 fim point source, while L906E-IRS 2 
is markedly fainter than L906E-IRS 1. IRS 1 is extremely 
red (particularly in comparison to L906E-IRS 2; Table 0, 
and is associated with extensive H2 filaments that proba- 
bly delineate at least two jets (discussed further in Section 



A5). If we also include the cluster of YSOs associated with 
L906E-IRS 3/IRS 4, then we probably have an evolutionary 
sequence, with L906E-IRS 1 representing the youngest and 
L906E-IRS 3/IRS 4 the most evolved region. 

The bright Spitzer sources and the clustering properties 
of the embedded stellar population are discussed in more 
detail in Paper II. 



3.2.2 The outflows 

In addition to the embedded young stars, Spitzer also de- 
tected the main DR21 outflow at 4.5 ^m. In high-mass 
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Figure 8. WFCAM H2 image of a portion of the DR21/W75 tile with, superimposed, contours of the 850 /xm emission. The SCUBA 
contours represent 0.125, 0.25, 0.5, 1, 2, 4, 8 and 16 Jy/beam (the beam size measures 14" at 850 /im). The extent of the SCUBA map 
is marked with dashed lines. 
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star forming regions like DR21/W75, distinguishing line- 
emission jet features from background diffuse dust and PAH 
emission regions is difficult, even in the relatively PAH-free 

4.5 ^tm band. With the broad IRAC bands one is not able 
to "continuum-subtract" the data. One must therefore use 
colour images as a diagnostic tool. We show such an image 
in Fig. □ 

As noted bv | Marston et al.l (|^004*) and more recently by 
ISmith et al.) l)200(ih . the main DR2f outflow is detected by 
Spitzer. Morphologically, the fl ow is very similar in images 
at 2.12 /xm and 4.5 /^m. Indeed. ISmith et al.l ll20 06l) overlay 
contours of the H2 l-OS(l) emission o f iDavis fc Sm 'ithl dlQodi l 
onto a 4.5 /im image and find an almost exact correlation be- 
tween the emission-line features in the two bands. Note also 
that the flow is consistently the same colour (pink) in Fig. 
Although the IRAC filter at 4.5 /im encompasses molecular 
and atomic emission lines from CO (v=l-0 at 4.45-4.95 pim) 
and HI Bra (at 4.052 /im) , the unchanging colour along the 
DR21 outfiow in Fig. |7|suggests that the bulk of the 4.5 /xm 
emission is in fact from the H2 0-0S(9) line (at 4.694 /xm). 
This conclusion is supported by ISO spectroscopy of por- 
tions of the DR2 1 outflow (Smith. Eisloffcl & Davis 1998; 
ISmith et al.ll200d) . 

In shock models, H2 emission in outflows is predicted to 
be particularly strong in the 4.5 lira band, 5-14 times higher 
than i n the 3.6 ^im band CS mith fc Rosenl2005^ . ISmith et all 
J2006I) find that tn outflows ~50% of the flux in the IRAC 

3.6 Atm, 4.5 Mm and 5.8 /xm ban ds is due to H2 line emission. 
ISmith. Eisloffel fc Davi^ lll99dl propose excitation of the H2 
emission in multiple C-type bow shocks. ISmith et all ll2006l) 
also note a conspicuous lack of CO, Bra and PAH emission 
in the ISO spectra. 

Most of the "red" H2 features evident in our WFCAM 
images (Figs. |5| and are detected by Spitzer, albeit at 
lower spatial resolution. In Fig. |7|one or two features along 
the ridge appear somewhat "redder" than other jet features, 
probably because of increased extinction (note in particular 
the knot ~1' north of DR21-IRS 4). This may also explain 
the brown colouration at the end of the western DR21 flow 
lobe, although Smith ct al. ( 2006) - who also note this subtle 
change in colour - mention that a reduced pre-shock velocity 
or increased pre-shock density may have the same effect. 

The Spitzer 5.8 /xm band covers emission from the 
strong [Fell] line at 5.34 /^m. However, upon examining the 
5.8 fj,ra images (not shown) we found it impossible to distin- 
guish line emission features associated with the outflow from 
the more intense, fllamentary 6.2 /xm PAH emission that per- 
vades the 5.8 /xm data (the 3.6 /xm and 8.0 fim bands likewise 
include bright PAH emission at 3.3 /^m and 7.7 /^m). Suffice 
to say that no bright [Fell] component was found in the main 
DR 21 outflow. 

3.3 Comparison with SCUBA data 

In Fig.|H|we overlay contours of 850 ^m emission on to our H2 
image. The 850 /xm scan map covers a region approximately 
12'x29' in size, encompassing both DR21, W75N and the 
region in between. A comparison of SCUBA and Spitzer 
8.0 ura. data is also made in Fig. |1| Parameters derived for 
the major cores in our map from Gaussian fltting of our flux- 
calibrated 850 /xm data are given in Table |5| The integrated 
flux is the sky-subtracted value within the elliptical aperture 



defined in the table. The relationship (if any) to 8.0 /im point 
sources is also noted, although we find measurable offsets 
between the SCUBA peaks and mid-IR sources in almost 
all cases. This is probably due to the fact that the SCUBA 
cores, and in some cases 8.0 /xm peaks, encompass multiple 
sources. Clearly, care needs to be taken when comparing 
source fluxes at mid-IR and far-IR wavelengths. 

The DR21 region comprises three groups of molecular 
cores, the DR21 Hll region itself, DR21(0H) and, a further 
3' to the north, a chain o f cores labeled FIR 1, 2 and 3 by 
IChandler. Gear fc ChinJ (^933) • New molecular cores are 
labeled FIR 4 to FIR 12, while sub-components within the 
DR21 and DR21(0H) regions are identified based on their 
relative locations. The brightest 850 /im peak, labeled DR21 
in Fig. |H| (and Fig. 2J , has no mid-IR c ounterpart, though 
it coincides with the Hll regions A-B-C iCveanowski et all 
.2003 ). The 850 /im peak towards DR21(0H) coincides pre- 
cisely with the 3.1mm con tinuum source MMl observed by 
iLiechti fc WalmslevI ^1997^ . 

As already noted, ERO 1 is detected at 850 /im. How- 
ever, neither ERO 2 nor ERO 3 coincide precisely with an 
850 /^m peak; ERO 2 in particular is offset ~15" from the 
nearest SCUBA core, DR21(OH)S. At 8.0 /xm ERO 4 com- 
prises three 8 ^.m sources. Close examination of Fig. |31 sug- 
gests that at least two of these emi ssion features - ERO 4c 
and ERO 4s in the nomenclature oflMar ston et al.l ^2004^ - 
are actually the illuminated bicones of scattered light on ei- 
ther side of a band of obscuration, possibly associated with 
an edge-on disk. A K-band point source sits within the apex 
of the brighter, northeastern cone, while a second, much red- 
der source occupies the southern cone. Notably, this second 
source is closest to the faint SCUBA core FIR 6, although 
again these features are not precisely coincident. 

There are additional cores detected at 850 /im that are 
not associated with bright 8.0 /im sources (Table HI . These 
may harbour very young sources or be pre-stellar cores. 
North of DR21(OH) and ERO 3, it is interesting to note 
that the 850 /im emission FIR 1/FIR 2/FIR 3 roughly tracks 
the chain of reddened stars rather than the dark patch where 
very few stars are seen (even with Spitzer). Even so, it seems 
clear that the 8.0 /im peaks and SCUBA cores are associated 
with different populations of young objects. The SCUBA 
cores may eventually evolve into compact clusters of young 
stars, much like those seen with Spitzer and WFCAM (the 
"IRS" sources), in which a dominant 8 /tm peak, probably 
an early B-type star (Paper II), is surrounded by a number 
of fainter, lower-mass sources. 



3.4 The H2 2.122 /im outflows and their driving 
sources 

For the sake of brevity we consign our description of the 
many H2 jets and outfiows in DR21/W75 to the appendix. 
In Figs. to El we show continuum-subtracted H2 im- 
ages of various regions along and adjacent to the DR21/W75 
molecular ridge. For simplicity we split the DR21/W75 re- 
gion into five areas, which we label A-E (see Fig[TJ. We 
then assign numbers to individual H2 features in each re- 
gion. Knots that appear to be part of the same jet are given 
consecutive numbers. For example, knots A 1-1, A 1-2 and 
A 1-3 are possibly part of the same jet in area A, while A 2- 
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Table 1. Catalogue of bright 8 fim sources in DR21/W75 



Rcfl 


(2000.0) 


Dcc^ 

t z U U U . u _) 


[3.6] 


[4.5] 


[5.8] 


[8.0] 




Name 


Descriptive comment 




20:39 


16.72 


42: 16 : 09. 1 


8.12 


7.06 


5.49 


4.82 


1.1 


EROl 


__ 

K-band/8^tm point source ("a" in Table^H 


2 


20:39 


02.86 


42:22:00.2 


10.48 


9.86 


7.29 


5.58 


3.2 


ER02 


Diffuse 8 /.tni extending ~5" W of K-band star 


3 


20:39:00.47 


42:24:36.6 


9.74 


8.00 


6.87 


5.93 


1.5 


ER03 


K-.band/8^im point source 


4 


20:39 


08.33 


42: 15:44.9 


11.25 


10.48 


9.32 


7.80 


1.2 


ER04c 


K-band point source in NB 8.0 ^tm arc 




20:39:07.28 


42' 15 ■ 34 8 


12.06 


11.38 


8.65 


6.90 


3.0 


ER04E, 


K-band point source in SW 8 /J.m arc 


6 


20:39:00.05 


42:19:36.6 


_ 


_ 


_ 


_ 


_ 


DR21~IRS1 


No distinct 8 fim peak 


7 


20:38 


56.91 


42:18:58.8 


_ 


_ 


_ 


_ 


_ 


DR21-IRS2 


Diffuse 8 Mm peak ^3"N of K-band star 


S 


20:38 


59.83 


42:18:58.1 


10.74 


10.18 


9.44 


8.53 


-0.3 


DR21-IRS3 


Faint K-band/8 ^m point source 


9 


20:38 


59.32 


42:20:16.1 


9.82 


9.14 


8.02 


6.33 


1.2 


DR21-IRS4 


Faint K-band/8 ^m source with extended emission 


10 


20:38 


55.85* 


42:21:54.1* 


12.06 


11.37 


8.32 


6.47 


4.1 


DR21-IRS5 


Diffuse 8 fim peak; no K-band point source counterpart 


11 


20:38 


57.19 


42:22:41.1 


7.40 


6.62 


5.35 


4.13 


1.0 


DR21-IRS6 


Bright K/S/zm point source ( "e" in Tabloid 


12 


20:39:02.00 


42:24:59.2 


9.22 


6.92 


5.46 


4.54 


2.4 


DR21-IRS7 


Bright 8 /xra source with faint K-band counterpart 


13 


20:39 


03.72 


42:25:29.7 


7.96 


6.62 


5.17 


4.46 


1.3 


DR21-IRS8 


Bright K./8 fim source 


14 


20:39:09.49 


42:24:26.8 


9.32 


8.51 


7.86 


7.03 


-0.2 


DR21-IRS9 


K/8 fim point source 


15 


20:38:46.37 


42:24:39.7 


7.63 


6.91 


5. 86 


5.25 


-0.0 


DR21-IRS10 


Slightly extended K/8/^m point source 


16 


20:38 


36.76 


42:37:29.4 












W75N-IRS1 


SE source of three 8 fim peaks 


17 


20:38 


35.37 


42:37:13.4 












W75N-IRS2 


Nebulous source at K and 8 fim. 


IS 


20:38 


38.82 


42:37:17.6 


8.96 


8.30 


7.65 


6.32 


0.2 


W75N-IRS3 


K-band/8 Mm point source 


19 


20:38 


37.73 


42:37:59.4 












W75N-IRS4 


Diffuse 8 fim emission enveloping K-band point source 


20 


20:38 


37.29 


42:39:33.0 


10.08 


9.62 


7.21 


5.44 


2.9 


W75N-IRS5 


Bright K/S^m point source 


21 


20:38 


40.16 


42:39:53.5 












W75N-IRS6 


K-band point source in shell of 8 ^J.m emission 


22 


20:38 


08.27 


42:38:36.3 


7.09 


6.26 


4.80 


4.26 


0.6 


W75N-IRS7 


K-band/8 Mm point source ( "m" in Tablc^ 


23 


20:38 


07.10 


42:38:52.1 




7.38 




5.39 




W75N-IRS8^ 


K-band/8 Mm point source 


24 


20:37 


55.40 


42:40:46.9 




8.27 




5.76 




W75N-IRS9'^ 


K-band/8 Mm point source 


25 


20:37 


47.48 


42:38:37.0 




9.71 




8.00 




W75N-IRS10^ 


K-band/8 Mm point source 


26 


20:37 


47.25 


42:38:42.0 




10.04 




6.95 




W75N-IRS11-^ 


K-band/8 Mm point source 


27 


20:37 


47.35 


42:40:53.0 




7.69 




7.19 




W75N-IRS12^ 


K-band/8 Mm point source 


28 


20:37 


57.77 


42:38:54.0 




11.06 




6.78 




W75N-IRS13'^ 


K-band/8 Mm point source 


29 


20:38 


38.69 


42:35:34.2 


11.69 


11.31 


10.21 


8.79 


0.6 


W75N-IRS14 


K-band point source in shell of 8 Mm emission 


30 


20:37 


17.67 


42:16:37.4 


10.62 


10.03 


6.83 


5.04 


4.1 


L906E-IRS1 


K/S Mm point source ( "n" in Tablc^l 


31 


20:37:26.35 


42:15:51.9 


8.58 


8.00 


7.44 


7.02 


-1.0 


L906E-IRS2 


K/8 Mm point source in shell of 8 Mm emission 


32 


20:37:27.73 


42:14:13.1 


8.72 


7.78 


6.53 


4.93 


1.6 


L906E-IRS3 


Bright K/8 Mm point source 


33 


20:37:30.50* 


42:13:59.2* 


11.01 


9.09 


7.67 


6.67 


2.1 


L906E-IRS4 


Bright 8 Mm point source; no K counterpart 


34 


20:37:31.91 


42:13:01.5 


11.64 


11.03 


8.39 


6.59 


3.4 


L906E-IRS5 


Bright K/8 Mm point source 


35 


20:36 


57.90 


42:11:29.7 




10.46 




8.92 




L906E-IRS6^ 


Nebulous K-band source(s) 


36 


20:39 


01.27* 


42:19:54.1* 


10.21 


7.52 


5.53 


3.98 


4.2 


DR21-D 


Very bright 8 Mm point source; no K counterpart 


37 


20:39 


01.01* 


42:22:50.2* 


11.80 


10.74 


8.13 


6.27 


3.8 


DR21(OH) 


8 Mm point source; no K counterpart 



Reference number used in the colour-colour diagrams in Paper II. Those arc not the "IRS" designations used here in the figures, 
^Coordinates of K-band point source counterparts to bright 8 /xm features in DR21/W75. Coordinates marked with a * have no K-band source, so the 8 /xm peak 
position is listed, 

^These targets were not observed in all IRAC bands. 



1, although still in area A, is probably part of a different 
flow. 

In Table El where possible, jet lengths, position an- 
gles and jet opening angles are estimated, although many 
emission-line features are single knots with no obvious cen- 
tral engine. We do not discuss the main DR21 outflow, 
since this H2 flo w has been studied in considerable detail in 
past papers (e.g. [G arden et al."l991a|; iDavis fc Smithlll996l : 
ISmith. Eisloffel fc Daviail998i) . In all we identify about 50 
separate outflows. 

Although in previous sections we have identified and 
labeled the bright K-band/8. /xm sources, these are not 
necessarily the outflow sources. We have therefore used 
the Spitzer IRAC images to identify young stellar objects 
(YSOs) based on their spectral index. As mentioned earlier, 
Q is derived from linear fits to the observed flux density in 
each of the four IRAC bands. Note that the more nebu- 
lous sources, like DR21-IRS 1 and W75N-IRS 1, were not 
retrieved by our analysis technique (described in Paper II) 
and therefore have not been flagged as YSOs. 

Traditionally, YSOs have been classified as Class 0, I, II 
or Class III based on their bolometric temperature, infrared 
spectral slope (measured between 2 fim and 20 n m) and the 
associated colour excess due to a dusty envelope llLadall987l : 
[Andre. Ward- Thompson fc Barsonv .1993.'l . Class SEDs 
consist entirely of grey-body emission from the cool cir- 
cumstellar envelope and are t hus detectable only at mid- t o 
far-infrared wavelengths (e.g. iNorieea-Crespo et al.l 12004) . 
Class I, flat-spectrum. Class II and Class III YSOs, de- 



tectable in the near-infrared, have photospheric SEDs with 
an infrared excess, the slope of which, a, is positive f or the 
younger (Class I) sources. More recentlv. iLada et alj ll200d) 
have applied this classification scheme to Spitzer photom- 
etry of YSOs in the low-mass star forming region IC 348. 
They distinguish stars with optically thick circumstellar 
disks from their "anemic" counterparts, sources with heav- 
ily depleted, optically-thin disks or disks with large central 
holes. Their thick-disk sources have SEDs measured from 
the IRAC bands of q > —1.80, while the anemic disk sources 
have shallower spectral slopes with —1.80 > a > —2.65. 

In principle, IRAC colour-colou r diagrams can also 
be used to distinguish YSQ class es (jMeeeath et aljEo04l : 
I Allen et al.ll2004l:lLada et al.l200d) . In the four IRAC bands 
centred at 3.6 /im, 4.5 /im, 5.8 /xm and 8.0 /xm. Class and 
Class I protostars have roughly the same colours; [3.6]- 
[4.5]>0.7 and [5.8]-[8.0]>1.0, indicative of mass accretion. 
Class II sources (T Tauri stars) overlap this region, roughly 
occupying a box with 0.0<[3.6]-[4.5]<0.8 and 0.3<[5.8]- 
[8.0]<1.1. Reddened Class II YSOs have somewhat higher 
[3.6]-[4.5] colours, although slightly lower [5.8]-[8.0] colours. 
However, as we discuss in Paper II, there is considerable 
overlap in colour-colour spa ce between the various YSO 
classes. FlVIegeath et al.l(|2004^ note that from Spitzer colours 
alone it is possible to mis-identify reddened Class II sources 
or Class II sources with edge-on disks as Class I sources. 
Also, strong silicate absorption at 9.7 /xm and water-ice ab- 
sorption features at 3 fim and (to a lesser extent) 6 fim may 
modify the Spitzer colours of massive young stars, increas- 
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Table 2. Characteristics of the brightest 850 /im cores in DR21/W75 



Source 


RA 


Dec 


Flux" 


Sizc^ 


P.A.^ 


8 /jm counterpart 




(2000.0) 


(2000.0) 


(Jy) 


(arcsec) 


(degs) 




DR21^ 


20 39 00.91 


42 19 35.9 


54.5 


27 xl6 


13 


No mid-IR counterpart 




90 "^Q no 7n 
oy uu- 1 u 


A 9 1 S 1 1 ^ 

^Z ±0 L L.O 


u.uo 


U.^ X o.u 


-oo 


No mid-IR counterpart 


riR 91 QTT 


90 "^Q fif^ QQ 


AO 1 8 1 Q A 
4Z lo iy.4 


'\ 1 
0. 1 


9Q V 91^ 

zy A zo 


-oy 


No mid-IR counterpart 


UJXZi VV i 


on 'iQ f^n 
ZU v5o 04. ou 


/I o 1 Q 1/1 n 
4Z ly 14. u 


D.O 


QQ v 1 Q A 
Oo X lo.4 


1 Q 

lo 


~ou VV IN VV 01 urtzi-irxo z, wesTiern uixzi now 


FIT? 91 '\A/9 
LJiXZl VV Z 


on 'iS '^(^ 


A9 1 s f^n n 
4Z lo ou.u 


n 9Q 


1 c: c; V Q 7 
lO.O X y . i 


9S 

-Zo 


No mid-IR counterpart 
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DiflFuse 8 /j.m emission; eastern DR21 outflow lobe 
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No mid-IR counterpart 
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ZU Oo oy.oz 
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o.D 


A9 V 1 9 
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1 (\ 
ID 


No obvious mid-IR counterpart 


FIR 1 


20 39 00.20 


42 24 37.0 


4.4 


29.8x8.3 


-85 


ERO 3; ~2.0" W of elongated 8 /^m peak 


FIR 2 


20 39 02.16 


42 25 01.3 


9.2 


22.5x11.9 
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~2.5" NE of elongated 8 /^m pea 
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No mid-IR, counterpart 


r irt 


on QQ 1 Q1 
ZU oy iO-oi 


A9 1 1 1 Q 
4Z ID 11.0 


'\ 9 
O.Z 


1 c; 1 V 1 9 1 
lO. 1 X IZ. 1 


- / 


ERO Ij coincides with pair of 8.0 /^iui sources 
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l.U 


Qfl v 1 Q 
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1 
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Hiixyj 4, VV o VV OI iv-DanQ/ o.u /im source HiIxkj 


7 
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ZU Oo 4D.04 


A O OA Q7 1 
4z Z4 / . 1 


n /in 

U.4U 
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( .O X D.O 
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-ZO 
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on Qo nn nQ 
ZU Oy UU. Uy 


A O OI QO n 
4z Zi OZ.U 


1 7 


OO V on 

zZ X ZU 




No mid-IR counterpart 


FIR 9 


on o o A r\ 1 r\ 

zU oo 49.10 


4z z7 40. o 


0.63 


ic n>^io o 

15. Ox 12. o 


-36 


No mid-IR counterpart 


FIR 10 


20 38 45.50 


42 28 10.3 


2.2 


52x14.8 


-22 


No mid-IR counterpart 


FIR 11 


20 38 43.76 


42 29 03.5 


0.62 


19.3x15.5 


-48 


No mid-IR counterpart 


FIR 12 


20 38 15.40 


42 32 44.8 


0.62 


16.7x15.2 


-59 


No mid-IR counterpart 


W75N9 


20 38 36.41 


42 37 34.2 


26.4 


11.8 x8.8 


10 


Coincides with elongated 8 ^m peak, to within ^^1". 


W75N-E 


20 38 47.28 


42 38 02.0 


1.1 


22x6.3 


-83 


No mid-IR counterpart 


W75N-N 


20 38 33.18 


42 39 45.7 


3.8 


31x27 


-22 


No mid-IR counterpart, though near outflow sources 


W75N-W 


20:38:31.02 


42:37:47.6 


2.3 


39x10 


81 


No mid-IR counterpart 


W75N-SW 


20 38 31.02 


42 36 27.1 


4.0 


43x33 


-85 


No mid-IR counterpart 


W75N-W1 


20:38:19.47 


42:36:52.5 


0.2 


11x7.6 


-71 


~5" S of * /im peak 


W75N-W2 


20:38:11.82 


42:37:34.6 


0.6 


21x3 


-20 


No mid-IR counterpart 


W75N-W3 


20:38:10.65 


42:38 04.6 


-l.O'^ 


9.1x8.7 




~40" SE of W75N-IRS 7 



"Integrated flux at 850/im. 

''Deconvolved size of the core and position angle of the long axis (me asured east of north). 
'^This primary peak is midway between Hll regions DR21 A,B and C ICvganowski et al.ll20oj) . 

''Probably includes CO J=3-2 line emission from the outflow; the SCUBA filter is centred at 353 GHz (passband ~30 GHz), while the 
rest frequenc y of CO 3-2 is 346 GHz. 

'^Resolved bylMangijm^^ootei^^^^m^J_^W2j) into two cores at 2.7 mm. 
.''Resolved bv lChandler. Gear fc Chinil ^1993a^ into two cores at 1.3 mm, Nl & N2. 

^Sub-mm core coinc ides with Hll region W75N{B), which is further resolved into three UC Hn regions VLAl, VLA2, & VLA 3 
iTorrelles et al.lll997^ . 

''Flux uncertain because this source is near the edge of the SCUBA map. 



in£_[3.6l-[4.5l though reducing [5.8]-[8.0] llGibb et alj|2000l: 
Ivan D ishoeck 2004). So although Spitzer remains a power- 
ful tool for distinguishing YSOs from background and/or 
foreground field stars, particularly in high-extinction re- 
gions where near-IR data are of limited use, here we only 
identify sources with reddened photospheres, and distin- 
guish protostars (a > 0.25) from pre-main-sequence objects 
(0.25 > a > -1.60). 

In Figs. lAll to IA5I we plot the positions of all of the 
YSOs that were extracted from the Spitzer images in all four 
IRAC bands, as reported in Paper II. We use circles to mark 
the protostars and triangles to indicate the remaining red- 



dened photospheres. Approximately 1600 sources were re- 
trieved across the entire DR21/W75 field of which 165 were 
protostars and an additional 482 were pre-main-sequence 
objects (we have only identified point sources with photo- 
metric errors of less than 0.05 mag at 8.0 fim). If we consider 
that some of the protostars will actually be Class II YSOs 
viewed edge-on, then the YSO population is reasonably con- 
sistent with an order-of-magnitude increase in the duration 
of the Class II phase over the Class O/I protostellar phase. 
The large-scale distribution and clustering of the YSOs in 
DR21/W75 is discussed in Paper II. 

The K-band magnitudes, mid-IR magnitudes and spec- 



Observations of jets in DR21/W75 



Table 3. H2 knot or outflow parameters (excluding the main DR21 and W75N flows) 



Jet 


RA" 
(2000.0) 


Dec" 
(2000.0) 


li- 

(arcsec) 


P.A.'^ 
(degs) 


Qd 

(degs) 


1 — 0S(1) 
(lO-i^W m-2) 


Descriptive note 


A 1 1 ■ A 1 Q 


on- QO- 1 
zu.oy. ly.z 


4Z. 14.0v5 


71 
1 1 


1 flO 

loy 


Q in 

0— lU 


on 
zu 


Oollimatea jet 


A 1 ■ A /I 


on- 'io- 1 7 n 
zu.oy. 1 ( .U 


/I O- 1 fl- 1 /l 
4Z. 10. 14 


07 
Z i 


^ ,A^ 
~40 


fin 
OU 


onn 

ZUU 


Compact, bipolar jet from ERO 1 


A Q 1 
A o-l 


zu.vSy.ui.y 


/I O- 1 Q ■ 1 1 
4Z. lO. 11 








QQ 
00 


Arc of H2 emission 
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/V 4-1, A. 4-z 
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zu.oo.Oo.y 
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4o 


OA 
Z4 
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zu 


Possible faint jet and bow 
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/V 0-1. Jt. O-O 
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ZU.OO.0D.4 


/I O- 1 fi- 1 fl 
4Z. lo. 10 


00 


00 

zy 


E^ Qn 
0— oU 


f^fl 
Oo 


Collimated jet 


A 1 ■ A fl 


zu. 00.41.0 


/I O- 1 o-ni 
4Z. ly .ui 








fin 
OU 


Possible bow shock 


A 7 1 
/A (-1 


ZU. 00. 40. 1 


/I O- 1 O- 1 /l 
4Z. ly . 14 








40 


Bright, extended, knot 


A C 1 
/V O-l 


on- Qfi-/i 7 
zu.00.4 i .0 


/I o- 1 o- 00 
4Z. ly. zz 








Ofl 
ZD 


Bright, diffuse knot 
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A 

4y 
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00 


Possible collimated jet 
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12 


A7 
4 I 


4U 


00 


Collimated jet 


Jj 1-1.13 1-0 


ZU.Oo.Oo. 1 


4Z.ZU.UO 


^^40 


111 
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A^ 
40 


AQ 
40 


Jet and bow shock 


ID Z-1 


on- Qfi- f^E^ n 
zu. 00. 00. U 


/I o-on - 








^ 1 n 
< lU 


Faint, elongated knot 
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zu. 00. oy.y 


/I O-Ofl- 
4Z. zu. 01 








1 K 
10 


Two or more fingers of emission 


B 4-1 


20:38:58.3 


42:21:09 


12 


175 


— 


39 


Extended bow or "bullet" 


B 5-1 


20:39:04.5 


42:22:44 
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Possible extension of B6-1;B 6-3 
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Bright, knotty jet 
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Curving jet 
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42:24:58 
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Possible faint/diffuse jet 


B 12-1:B 12-2 


20:38:45.9 
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4z.oy.uy 








E^n 

OU 
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0-1.15 0-0 
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42:37:12 
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44 
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Collimated jet 
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Faint chain of knots 
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Two faint knots (possible jet) 
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4z.oy . iz 








---in 
<.1U 
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---in 
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D 10-1 


20:38:03.2 
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_ 


_ 




<10 


Faint, compact knot 


D 11-1 


20:37:52.3 


42:40:50 








<10 


Faint knots and filaments 


D 12-1:D 12-3 


20:37:44.7 


42:40:27 


~50 


~110 




<10 


Faint group of knots (possible jet) 


E 1-1 


20:36:59.8 


42:11:19 


24 


118 


1 ! 1 


<10 


Faint knotty jet 


E 2-l:E 2-3 


20:37:02.3 


42:11:51 


21 


~139 


18 


23 


Possible knotty jet 


E 3-l:E 3-4 


20:37:00.1 


42:13:14 


~50 


~92 


40 


80 


Curving knotty jet 


E 4-1 


20:37:21.0 


42:16:03 


~40 


124 


30 


36 


Knotty jet 


E 5-1:5-6 


20:37:28.9 


42:16:12 


130 


~115 


12 


64 


Extended, knotty jet 


E 6-1:6-2 


20:37:19.0 


42:11:39 


40 


95 




<10 


Two faint knots 


E 7-1 


20:37:14.5 


42:13:01 








<10 


Faint bow shock 


E 8-1 


20:37:05.2 


42:13:14 








<10 


Faint knot (part of E 3-l:E 3-4?) 


E 9-1 


20:37:32.0 


42:14:27 


~25 


~175 




<10 


Faint north-south filament 



"Coordinates of the first H2 feature in each jet, A 1-1, A 2-1, B 1-1, etc. (sec Figs. IXTl to lA5i . 
''Length of the H2 flow, when multiple components are identified. At 3 kpc, 1" = 0.0145 pc = 3000 AU. 
'^Jet position angle (E of N), again, when multiple components are identified. 
''Flow opening angle (see text for details). 

"^Integrated H2 l-OS(l) flux, uncorrected for extinction, accurate to 20-30%. 
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Table 4. Photometry of possible H2 jet sources (excluding the main DR21 and W75N flows) 



Jet 


Source 


RA 


Dec 


K 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


a 


Other 




Ref " 


(2000.0) 


(2000.0) 














name 


A 2-l:A 2-4 


a 


20:39:16.73 


42:16:09.2 


10.98 


8.12 


7.06 


5.50 


4.82 


1.1 


ERO 1 


A 3-1 


b 


20:39:01.48 


42:18:09.7 


15.37 


14.55 


12.88 


11.88 


10.53 


1.7 




A 4-l:A 4-2 


c 


20:39:01.10 


42:18:59.0 




15.14 


13.17 


11.95 


10.27 


2.6 




A 5-l:A 5-3 


d 


20:38:57.35 


42:18:32.9 


14.67 


13.85 


13.11 


11.29 


9.60 


2.2 




B 7-l:B 7-3 


e 


20:38:57.20 


42:22:41.1 


8.34 


7.40 


6.62 


5.53 


4.13 


1.0 


DR21-IRS 6 


B 6-l:B 6-3 


f 


20:39:01.90 


42:23:02.3 




14.82 


13.57 


11.87 


10.26 


2.5 






g 


20:38:59.97 


42:23:07.7 


19.6 


15.47 


13.77 


12.30 


10.70 


2.6 




B 9-1 


h 


20:39:04.06 


42:26:07.8 


15.22 


13.22 


11.36 


10.46 


9.04 


1.8 




C 8-l:B 8-5 


j 


20:38:33.30 


42:39:38.5 


15.92 


12.51 


10.63 


9.26 


8.06 


2.2 






k 


20:38:32.88 


42:39:35.7 




13.02 


11.09 


10.04 


9.16 


1.5 




C 9-l:B 9-2 


1 


20:38:32.32 


42:39:45.2 


14.18 


11.19 


9.70 


8.76 


7.80 


1.0 




D 3-l:D 3-2 


m 


20:38:08.27 


42:38:36.4 




7.09 


6.26 


4.80 


4.26 


0.6 


W75N-IRS 7 


E 4-1 


n 


20:37:17.66 


42:16:37.4 




10.62 


10.03 


6.83 


5.04 


4.1 


L906E-1RS 1 



"Possible source of the jet identified in column 1 (see appendix for details). 



tral indices of the candidate outflow sources are given in 
Table H Most are undetected at the shorter near-IR wave- 
lengths; three were undetected even at K (fainter than ~19). 
Two of the brighter targets were saturated in our K-band 
data. The sources are assigned a letter in Table|l]and in each 
of the Figures in the appendix (note that some sources have 
both an IRS designation and an outflow source reference 
letter) . 



4 DISCUSSION 

4.1 Properties of the H2 jet sources 

In Fig. |5|we plot mid-IR spectral index, a, against source 
magnitude at 8.0 /im for all point sources in the DR21/W75 
field detected in all four IRAC bands. The plot suggests 
a mild correlation between source magnitude and spectral 
index. In other words, many of the bright 8.0 ^m sources 
appear to be embedded. The brighter 8.0 sources are 
therefore probably luminous young stars. Accurate luminos- 
ity estimates for individual sources are difficult to measure 
without photometry at longer wavelengths: note that the 
brighte st sources were satur ated in the MIPS images of this 
region ijMarston et alJl2004 ). and that IRAS data lack suffi- 
cient spatial resolution to distinguish individual stars. Even 
so, it seems safe to assume that the intermediate and high- 
mass sources correspond with the reddest 8.0 fim sources in 
Table El 

We find little evidence for H2 jets associated with the 
8.0 fim "IRS" sources (see the Appendix for a detailed dis- 
cussion), and only one of the four EROs has a definite H2 
jet. Most of the candidate jet sources in Table |1| are rel- 
atively faint at 8.0 ^m, while ~80% of the outfiows listed 
in Table El have no obvious source identified at 8.0 fim. It 
thus seems clear that, with the exception of the massive 
DR21 and W75N outfiows themselves, the H2 jets are not 
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13 01 

Figure 9. Plot of spectral index, a, against source magnitude at 
8.0 /^m. The sources listed in Table 111 are marked with squares, 
while the remaining IRAC sources are marked with crosses. 

driven by the massive protostars. Instead, low- or at most 
intermediate-mass YSOs seem to power each flow. 

The candidate outflow sources are, however, very red 
(in comparison to the criteria used to select protostars from 
the Spitzer observations). Many will probably be Class or 
early Class I sources. Some may be flat-spectrum or Class 
II YSOs viewed edge-on and therefore through their circum- 
stellar disks, although if this were the case one would expect 
to see more Class II jet sources (with negative values of a) 
viewed pole-on. It seems likely that many of the H2 flows are 
powered by low-mass "protostars" that are too faint and/or 
too heavily embedded to be detected even by Spitzer. 

Candidate central engines for the massive DR21 and 
W75N outflows remain elusive. Accurate Spitzer photome- 
try of the IRS sources in each region is difficult, since at 
near- and mid-IR wavelengths these sources are often ex- 
tended and bathed in diffuse emission. The brighter sources 
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are also saturated at the longer wavelengths. In DR21, 
ISmith et all (|200l) note that UCHll region DR21(D) is 
detected in the Spitzer bands and has an extremely red 
mid-IR SED. Even so, it is not well placed mid- way be- 
tween the extensive H2 flow lobes. Also, it does not coincide 
with the brightest dust-continuum peak in our 850 /xm data. 
DR21(0H), on the other hand, is associated with a SCUBA 
core; it is also detected at mid-IR wavelengths, and again 
is extremely red. Its outflo w is clearly tra c ed in methanol 
maser emission Jpiambcck fc Mentenlll99ol : iKoean fc sivshl 
Il998l. see also Fig. IHJ, though DR 21(OH) itself is as- 
sociated with multiple water ma sers iForster et alJ llQTSt 
iManeum. Wooten fc Mundvlll992h and therefore probably 
comprises multiple sources, any one of which could drive 
the outflow. 

The relationship between the infrared sources and mul- 
tiple UCHlI regions in W75N are discussed in detail by 
IShepherd. Kurtz fc TestH ^2004^. The mid-IR properties o f 
these sources are analvsed bv lPersi. Tapia fc SmithI i2006l) . 
From their analysis of ground-based and Spitzer mid-IR 
photometry they propose that W75N-IRS 2 is a B3 star, 
while W75N-IRS 3 and W75N-IRS 4 are intermediate- 
mass Class II YSQs. In Fig. El W7 4 N-IRS 1 is clearly 
elongated; here iPersi. Tapia fc SmithI ||200^ identify two 
sources, one coincident with the near-IR W75N-IRS 1 peak, 
and a second, much redder source coincident with the 
MMl dust continuum pea k (Shepherd. Testi fc Starkll200^ 
and VLA radio sources jTorrelles et al.lll997^ .' This sec- 
ondary mid-IR source is almost certainly associated with the 
W75N outflow source. Its position, RA(2000);20''38'"36.5" 
DEC(2000):42°37'33.6"is also notably coincident, to within 
an arcsecond, with the "W75N" 850 /im peak in Table. |5| 

Overall, W75N and DR21 are similar in many ways. 
Both regions harbour a massive core containing clusters of 
Hll regions, some of which are ultra-compact. Both contain 
water, OH and methanol masers. Both are also associated 
with bright, near- and mid-IR nebulosity and clusters of tens 
of embedded sources, many of which are B-type (5-10 M0) 
stars. Both are also surrounded by smaller, more compact 
clusters of just a few (or just one) early B-type stars plus as- 
sociated low-mass stars. DR21 and W75N themselves drive 
very massive molecular outflows, which are almost parallel 
on the sky (as is the flow from DR21(OH)). And flnally, nu- 
merous smaller-scale H2 jets and outflows are found in the 
vicinity of each massive star-forming region, these jets being 
driven by low or at most intermediate-mass young stars. 

4.2 General outflow characteristics: sizes, opening 
angles and orientations 

Estimated jet parameters for the more prominent outflows 
are given in Tabled At a distance of 3 kpc a parsec-scale jet 
should be >70" in length. At least five of the H2 outflows 
meet this requirement. For a modest flow velocity of 100 km 
s~^ a jet should reach a length in excess of a parsec in ~ 
10* yrs; thus, a Class source could easily be associated 
with a parsec-scale flow. Some of the long jets in Table |21 
like B 6-l:B 6-3 and C 8-l:C 8-5 do have candidate sources 
with SED slopes consistent with very red YSOs (Table [IJ. 

Precise jet opening angles (9) are rather difficult to 
meeisure, particularly at kilo-parsec distances. The emit- 
ting flanks of large-scale bow shocks are often much wider 



than the underlying jet (HH 212 is a spe ctacular example; 
IZinnecker. McCauehrean fc Ravnerlll996l) . while changes in 
flow direction over time will increase the apparent flow open- 
ing angle. Moreover, if one includes all features along the 
flow axis, or just the extent of the leading bow shock, very 
different values can be inferred. Consequently, for a few flows 
in Tablel^we list a range of values; the smaller is measured 
between the candidate outflow source and the knot furthest 
from the source, while the larger value represents the maxi- 
mum apparent flow opening angle that includes all emission 
features. If no YSO is identified, the angles are measured 
between the two most distant knots in the jet. Overall, the 
jets in the DR21 and W75N regions appear to have high 
degrees of coUimation, consistent with the apparent youth 
of the candidate central engines in Table |3] 

Flow position angles on the sky are also listed in Ta- 
ble 01 The two main DR21 and W75N outfiows have al- 
most the same position angle (~60°-70°). The DR21(0H) 
molecular outflow, which is not obv iously traced in H2, has 
a sim ilar "east-west" position angle iLai. Girart fc CrutcheJ 
l2003h . as does the prominent H2 jet B 6-1 to B 6-3. Each 
of these flows is orientated roughly orthogonal to the chain 
of 850 /im cores and the molecular ridge that runs north- 
south through DR21 towards W75N. If we simply count the 
number of flows in regions A and B (including the main 
DR21, DR21(OH) and W75N outflows) with position an- 
gles in the range 45°-135° and compare this to the num- 
ber of flows with angles in the range 0°-45° and 135°-180°, 
we find a ratio of 11 to 6. (For the whole region imaged, 
the ratio is 25:10.) Moreover, a Kolmogorov-Smirnov test 
yields a probability of only 53% that the distribution of 
angles in areas A and B is homogeneously distributed be- 
tween 0° and 180°. It is therefore possible that the fiow 
position angles are to some degree dictated by the large- 
scale cloud morphology. If the north-south elongation of the 
DR21/W75 ridge is produced by gravitational collapse along 
east-west orientated magnetic fields lines, then the random 
motions of individual star forming cores within the ridge 
are only moderately successful at scrambling the polar axes 
of the outflow sources. This process may be even less ef- 
fective with the very massive YSOs, DR21, W75N, ERO 1 
and DR21(0H), where the fiows clearly are orthogonal to 
the north -south molecular ridge that extends through DR21. 
Note that lVallee fc Fieg3 (jioOg) find that the magnetic field 
direction, inferred from 850 /im linear polarimetry observa- 
tions, is orientated roughly orthogonal to the DR21 filament, 
as expected. 

4.3 Spatial distribution of outflows and their 
relationship to YSO clustering 

In busy star forming regions like DR21 and W75N, interac- 
tions between neighbouring stars as they pass by each other 
may inhibit disk formation, or even strip stars of their cir- 
cumstellar disks. Since accretion drives outfiows, YSOs with- 
out actively accreting disks will obviously not produce H2 
jets. In Figs. IAll to lA5l one can see that clusters of YSOs do 
not generally coincide with an abundance of H2 features. Al- 
though these areas are not entirely devoid of H2 jets, H2 fea- 
tures do seem to be more abundant in other regions, which 
show a paucity of YSOs (as identified by our Spitzer anal- 
ysis), e.g. the region '--^l' northwest of the western DR21 
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Figure 10. Plot of the number of YSO neighbours within a 50- 
pixel (20") radius of protostars that are thought to drive out- 
flows (thick/blue line) and those that have no obvious H2 jets 
(thin/green line). 



outflow lobe (Figs. IXTl and and 3' -5' west of the busy 
L906E-IRS 3/fRS 4/IRS 5 cluster (Figs.EHJ. A stellar den- 
sity plot is shown in Paper II. 

To investigate this issue further, we consider whether 
outflow sources have fewer neighbours than YSOs with no 
obvious H2 jet. To do this we count the number of neigh- 
bours each YSO has within a radius of 50 pixels (20", or 
0.3 pc). YSO neighbours are classified as Spitzer-identified 
protostars (a > 0.25) or pre-main-sequence stars (0.25 > 
a > —1.60), since all of these sources are likely to be in 
the same star forming regions and therefore potentially in- 
volved in star-star interactions or disk stripping. However, 
because all candidate outfiow sources appear to possess red- 
dened photospheres (e.g. Table llj, we count neighbours 
only for the protostars. We have confined the analysis to 
a 8.0' X 10.7' (1200x1600 pixel) region centred near DR21- 
IRS 5, a 8.0'x6.7' region centred on W75N, and a 8.0'x6.7' 
region centred midway between L906E-IRS 3 and L906E- 
IRS 6. In other words, we focus only on regions where H2 
outflows are abundant (except for the region west of W75N, 
where we lack complete Spitzer data and therefore source 
classifications) . 

The resulting histograms are plotted in Fig. [TOl The 
mean number of neighbours for candidate outflow sources is 
1.4, while the mean for the remaining protostars is 1.8. Both 
values are, however, skewed by the large number of sources 
with no neighbours (the N=0 column in Fig. llOt . Many of 
these sources lie on the periphery of the north-south DR 21 
ridge, and around the edges of the W75N and L906E cores. 
If we exclude the sources with no neighbours from the statis- 
tics, then the outfiow sources have a mean of 2.0 neighbours, 
while the remaining YSOs have a slightly higher mean of 
2.9 neighbours. The latter is evident as a secondary peak 
in Fig. 1101 The slight increase in the number of neighbours 
associated with non-outfiow YSOs is an extremely tentative 
result, not only because of the modest number of sources 
considered, but also because of uncertainties in identifying 
outflow source candidates. 

At first sight it seems that very few YSOs in the region 
drive detectable H2 fiows. By assuming that individual H2 



knots are parts of extended fiows, we identify only 53 H2 
outflows in Table |21 while our analysis of the Spitzer pho- 
tometry reveals 165 protostars spread over the entire 0.8° 
field. Some of the Class O/I protostars will actually be pre- 
main-sequence Class II sources seen edge-on, and we do not 
include some of the very faintest knots and wisps of emission 
in Table 13 Moreover, (1) the distance to DR 21 will limit 
our ability to detect some of the "weaker" flows, (2) extinc- 
tion ( approaching 100-200 magnitudes at V in the densest 
cores: iMoore et alJll991bl: IChandler. Gear fc Chiii] Il993al: 
IPersi. Tapia fc Smithll2006ft will hide some H2 features, and 
(3) east-west orientated outflows will probably break out of 
the denser cloud regions and therefore not be particularly 
bright in H2 emission. Even so, we conclude that we have de- 
tected - within a factor of a few - roughly the same number 
of H2 outflows as embedded protostars in DR21/W75. 

4.4 Momentum injection and cloud support 

Supersonic turbulent energy in giant molecular clouds is 
expected to decay on timescales of a few million years 
iOstriker. Stone fc Gammid 1200 j) . However, clouds last 
longer than this, harbouring multiple epochs of star for- 
mation at any one time. Outflows may be responsible for 
the continued injection of turbulent energy in to molecular 
clouds, maintaining virial equilibrium and thereby limiting 
cloud collapse and the overall star-formation efficiency. 

The turbulent momentum (and energy) in a cloud is 
approximately equal to the product of the cloud mass and 
turbulent line width (and line width squared). The mass of 
the cores associated with DR21/W75 have be en estimated 
from modeling of dust continuum obse rvations jMoore et alJ 
ll991al:IChandler. Gear fc "Chin][i9933), molecula r line maps 
JPickel et alll978HWilson fc MauersbergeJl990l) and, more 
recently, as virial masses from extensive CS observations 
fShi rlev et al.l2 003'V Core masses for DR21, DR21(OH) and 
W75N of between 1 x 10^ Mq and 3 x lO'' M© are pre- 
dicted. If we adopt a turbulent velocity of ~3-5km s~^, 
estimated from sub-arminute-resolution maps in gas tracers 
that are optically thin and/o r probes of high-density cores , 
like CS, NH3 and C^^O (e.g.lWilson fc MauersbergeJll99d: 
IChandler et"ai]ll993bl: IVallee fc Fiegdl2006^ - observations 
that are likely to yield velocities relatively unaffected by out- 
fiows or large-scale cloud motions - then a total momentum 
in the range 3 x 10^ - 1.5 x 10^ Mq km s~^ and turbulent 
energy in the range 1.8 X lO"*" - 1.5 X 10**^ J is inferred for 
each region. 

The momentum and energy supplied by an outflow from 
a young star vary with time (e.g. lSmithl2000l) . Younger fiows 
(from Class sources) are more powerful than their Class I 
and Class II counterparts, and therefore provide a larger mo- 
mentum fiux, P, and energy fiux (or mechanical power), L. 
.Walaw cnder. Ballv & Rcipurth (2005) note that one could 
use the canonical values for P for Class and Cl ass I so- 
lar mass stars measured bv lBontemps et alj il996l) to esti- 
mate the momentum supplied over the accreting lifetime of 
a young star. If P(Class 0) ~ 6 x lO"'^ Mq km s"^ yr"^ 
and P(Class 1) ~ 4 x 10~^ Mq km yr~^, then assuming 
durations of 10* and 10^ years for the Class and Cleiss 
I stages, over its lifetime an outfiow from a low-mass pro- 
tostar should inject ~1.0 Mq km s~^ of momentum into 
its surroundings. The outfiow kinetic energy scales with the 



Observations of jets in DR21/W75 19 



velocity; since most of the molecular material in outflows 
is probably entrained in the relatively slow-moving wings 
of bow shocks (which explains why CO outflows are usu- 
ally an order of magnitude slower than Hlf jets), then for 
a molecular flow velocity of ~20km s~^, the kinetic energy 
supplied by a typical outflow would be ~ 4 x 10'^^ J . The 50 
or so flows from low-mass stars in DR21/W75 could there- 
fore provide of the order of ~50 M0 km s~^ of turbulent 
momentum and inject about 2x10"^'' J of kinetic energy into 
the cloud over their lifetimes. Both values are less, by two to 
three orders of magnitude, than what is apparently needed 
to maintain the observed turbulent momentum and energy 
in the region. Note also that the trend toward east- west ori- 
entated outflows noted earlier, and the small filling factor 
of each coUimated jet (with respect to the total volume of 
the cloud) limits the efficiency of momentum and energy 
transport to the global GMC environment. 

But we have not considered the high-mass sources. Al- 
though these are relatively few in number, their outflows 
are more massive and more energetic, by three to four or- 
ders of magnitude. The momentum in the DR21 and W75N 
outflows, measured from CO observations, is of the order of 
~ 2 X lO'^ M(7)km s~^ and 1 0^ - lO'^ Mp, km s ~ \ respec- 
tively (|Garden et al.1 Il991bt IPavis et al.1 ll998aL ISheoherdl 
120011) . This would be more than enough to provide turbu- 
lent support against cloud collapse, if this momentum was 
distributed throughout the cloud. But of course these two 
flows do not impinge on the whole cloud; both are fairly well 
coUimated and so probably efficiently "drill" through their 
molecular surroundings, rather than feed turbulence into the 
larger environment. Note that massive dense cores, traced 
at 850 /im, remain around DR21, DR21(OH) and W75N 
(Fig. |5J; they are not entirely dispersed by the powerful 
outfiows. 

Finally, rather than use "mean" values for the momen- 
tum and energy of the outflows, we can infer these directly 
from the integrated H2 line fluxes, since the outflow me- 
chanical power should be directly re lated to the ene r gy ra- 
diated in the entraining shock fronts. iDavis fc SmithI lll996ll 
point out that in strong-shock prompt entrainment scenar- 
ios, where the ambient gas is swept up in a momentum- 
conserving fashion to form the molecular "CO" outflow, 
the energy radiated in the wings of the entraining bow 
shock - the flux seen in H2 and other coolants - should 
be equivalent to the mechanical power in the swept-up out- 
flow. In regions A, B and C, the total H2 l-OS(l) line 
flux, -Fi_o(si)i measured from the H2 flows identifled in 
Table El (again not including the special case DR21 and 
W75N flows) is about 1 x lO"^'' W m"^. At a distance 
d ~3 kpc, and assuming that the l-OS(l) line represents 
10% of all H2 line radiation, then the total H2 luminosity 
win be Lh2 = 10 X I/i_o(si) x Anif = 1 x lO^'^ W. Since 
ther e will be other important molecular coolants; H2O, CO, 
etc. jBenedettini et alJ 1200(1 |g iannini. Nisini fc Lorenzettil 
|20o3), the energy radiated in the outflow bow shocks will 
actually be 2-3 times higher. The observed line fluxes are 
also not corrected for extinction. The observed jets proba- 
bly occupy regions of fairly modest extinction; an of 10 
would boost the luminosity by a further factor of 2.5. The 
total radiated power is therefore probably an order of magni- 
tude higher, ~ 10^* W. For a dynamical age of ~ 3 x 10^^ sec 
(1 pc/100 kms"^) the energy pumped into the cloud equates 



to ~ 3 x lO''-' J, which a gain seems to be insufficient to ac- 
count for the turbulent energy in the GMC. 

In conclusion, the current generation of low-mass YSO 
flows probably does not drive the observed cloud turbulence. 
However, many repeats of the current epoch of outflow activ- 
ity, possible during the few million year timescale for turbu- 
lent decay in the GMC, might provide sufficient momentum 
and energy. Notably, the 10"^ — 10* year outfiow dynamical 
timescale is a factor of 10-100 times less than the turbulent 
decay time scale (for a ~3 pc cloud and velocity of ~3 km 
s~^), as required. Imaging surveys of other high-mass star 
forming regions are required to establish whether the abun- 
dance of outfiows in DR21/W75 is common in other regions, 
or whether DR21/W75 is uniquely active. The former would 
support the occurance of repeat epochs of outflow activity 
in this and other regions. 



5 SUMMARY AND CONCLUSIONS 

We compare near-IR and mid-IR images of the high-mass 
star forming regions DR21 and W75, obtained with the 
wide-fleld camera, WFCAM on the U.K. Infrared tele- 
scope, and the IRAC camera on the Spitzer Space telescope. 
850 nm images of much of the region, obtained with the 
SCUBA bolometer array on the James Clerk Maxwell Tele- 
scope, are also considered. 

The H2 l-OS(l) images obtained with WFCAM reveal 
numerous line-emission features, most of which appear to be 
associated with coUimated flows. The Spitzer IRAC images 
were found to be of limited use in tracing these outflows, 
because of the intense PAH emission in the IRAC bands as- 
sociated with background nebulosity, and because of limited 
spatial resolution in these very busy star- forming regions. 
However, the IRAC data were extremely useful for detect- 
ing YSOs, particularly massive, embedded sources, and for 
identifying potential outflow central engines. The near- and 
mid-IR photometry and stellar population are discussed in 
Paper II. The conclusions reached in this paper (Paper I) 
are described below. 

(i) The well-known W75N and DR21 high-mass star 
forming regions are surrounded by smaller clusters that har- 
bour perhaps just one intermediate or high-mass star (the 
DR21-IRS and W75N-IRS sources). In the case of DR21, the 
SCUBA data reveal less evolved cores that are not gener- 
ally associated with the "IRS-type" clusters, although these 
cores may eventually evolve into such clusters. The fact that 
the SCUBA cores are already distributed along the DR21 
ridge suggests that the cores are dynamically distributed 
throughout the region before they evolve into intermediate 
or high-mass star forming clusters. As noted in Paper II, 
the more evolved YSO population traced by Spitzer is even 
more widely distributed. 

(ii) Although we flnd a large number of H2 emission-line 
features, many of which appear to be associated with col- 
limated jets, these are usually not associated with either 
the young SCUBA cores or the compact infrared clusters 
noted above. Instead, they seem to be driven by independent 
protostars, probably low-mass Class O/I sources distributed 
around the periphery of the dense 850 /im cores and north- 
south molecular ridge. 
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(iii) H2 features associated with at least 50 individual out- 
flows are identified. Many are knotty and well coUimated, 
much like jets in nearby low mass star forming regions. At 
least five are parsec-scale flows. The orientations of the flows 
on the sky show some degree of order, with two thirds of the 
flows being roughly perpendicular (within a 90° cone) to the 
north-south molecular ridge; we find only a 53% chance that 
the flows are randomly orientated. 

(iv) We find very tentative evidence that clustering may 
inhibit disk accretion and thus the production of extensive 
outflows. We certainly do not see an enhancement in the 
number of outflows around compact clusters or groups of 
embedded protostars. 

(v) Although the low-mass YSO flows arc abundant and 
widely distributed throughout the region, the observed out- 
flow activity needs to be repeated 10-100 times if the flows 
are to provide sufficient momentum and kinetic energy to 
account for the turbulent motions in the large-scale molecu- 
lar cloud. The current high-mass YSO flows associated with 
DR21 and W75N are sufficiently energetic. However, be- 
cause these massive flows are well coUimated, orientated or- 
thogonal to the clumpy GMC that runs through DR21/W75, 
and because they arc few in numbers, the transport of their 
momentum and energy to the ambient medium is probably 
relatively inefiicient. 

The mid-IR and far-IR observations very effectively 
trace the distribution of protostars and pre-stellar cores 
across the entire region. The narrow-band H2 images com- 
plement these data by showing regions that arc dynamically 
active, areas that may otherwise go unnoticed. Together, the 
UKIRT/WFCAM, Spitzer/IRAC and JCMT/SCUBA data 
yield a near-complete picture of star formation across the 
extensive region targeted here. Similar wide-field studies in 
the future at modest to high spatial resolutions - partic- 
ularly with the advent of SCUBA-2 - should continue to 
improve our understanding of both low and high-mass star 
formation on global scales. 
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APPENDIX A: DESCRIPTION OF THE H2 
JETS 

In Figs.lXTltolXsl we show continuum-subtracted H2 images 
of various regions along and adjacent to the DR21/W75 
molecular ridge. We have separated the data into five ar- 
eas for clarity, labeled A-E in Fig. Individual H2 knots, 
filaments and bow-shocks within each area have been num- 
bered. Features that are thought to be part of the same 
outflow have been given consecutive numbers. For example, 
knots A 1-1, A 1-2 and A 1-3 are possibly part of the same 
jet in area A, while A 2-1, although still in area A, is prob- 
ably part of a different flow. We also mark on each figure 
the locations of Spitzer sources with reddened photospheres; 
those with a > 0.25 are dubbed protostars, and those with 
0.25 > a > —1.65 we refer to as pre-main-sequence objects. 
Note that only sources detected in all four Spitzer IRAG 
bands were extracted from the data. 



Al DR21 south - region A 

The area just to the south of DR21 is shown in Fig. lAll The 
region includes the 850 /xm core and mid-IR peak ERG 1, 
and the possible edge-on disk system ERO 4 discussed ear- 
lier. ERG 1 (star "a" in Table 0J is identified as a reddened 
protostar in our analysis of the Spitzer IRAG data (it is 
marked with an open circle). Bright H2 emission is associ- 
ated with ERO 1, though not with ERO 4. Knots A 2-2, 
A 2-3 and A 2-4 outline a roughly conical structure with 
the apex facing ERO 1. This morphology could be indica- 
tive of a wide flow opening angle. The counter-fiow is traced 
by knot A 2-1. The fiow axis is clearly orientated roughly 
northeast-southwest. 

A rmmber of coUimated H2 jets are found scattered 
around the region. The jet described by A 1-1, A 1-2 and 
A 1-3 is almost certainly unrelated to ERO 1. This jet has 
no Spitzer-detected central engine (no open circle or triangle 
in Fig. lAH . Nor is it associated with an 850 /xm continuum 
peak. 

Jet A 5-1 to A 5-3 and jet A 10-1 both point back 
roughly towards the infrared cluster associated with DR21 
(DR21-IRS 1 in Fig. Inland Fig.|XlJ, and could therefore be 
driven by sources within this region. Alternatively, A 10-1 
could be associated with the lu minous, accreting i nfrared 
source and Hll region DR21 D (|s^ith_iiii]|2QQi), while 
knots A 5-1 and A 5-2 could be driven by the bright K-band 
source DR21-IRS 2 (FiglU or the reddened YSO "d" found 
midway between knots A 5-1 and A 5-3. 

A 4-1 and A 4-2 possibly represent a faint jet and bow 
shock, with YSO "c" as a tentative source candidate. No- 
tably, A 4-l:A 4-2 is parallel with A 5-l:A 5-3. 

A 9-1, A 9-2 and A 9-3 may delineate another coUimated 
fiow that is roughly parallel with the main DR21 outflow, 
while knot A 6-1 may be a bow shock driven by a source 
to the southeast. Protostellar outflow source candidates are 
not found in the Spitzer data, however. 
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Figure Al. Continuum-subtracted H2 l-OS(l) image of region "A", just to tlic south of DR21. Inset are contour-plots showing the 
prominent line-emission features in more detail. The micro-stepped data have been binned over 2x2 pixels to a scale of 0.4" to improve 
the signal-to-noise on the extended H2 features. Contour levels measure 0.4, 0.8, 1.6x10-1* W m ^ arcsec ^ (black) and 3.2, 6.4, 
12.8x10-1* W m ^ arcsec ^ (white). The open circles mark the locations of the protostars {a > 0.25); the triangles mark the positions 
of the pre-main-sequence objects (0.25 > a > —1.60). Those that are possible jet sources are given ref erence letters and are also 
marked in the contour plots. The filled squares mark the positions of methanol masers in the DR21 outflow iPlambeck fc Mentenlll990t 
iKogan fc Slvshlil998i) : the star indicates the location of the DR21 A-B-C Ho region (DR21 D is detected by Spitzer and is labeled). 



A 3-1 is a fairly bright though non-descript feature in 
the centre of the field. It is close to the Spitzer YSO star 
"b" , which is a protostar with a notably steep spectral index 
(Table EJ. 

Additional faint knots and features are identified 
around the region which may or may not be associated 
with collimated flows. In many cases, the driving source of 
each flow is likely found along the molecular ridge that runs 
north-south through DR21 (Fig.|SJl. 



A2 DR21 north - region B 

Area B (Fig.|X2ll is centred on DR21(0H) and the jet labeled 
B 6-1 to B 6-3. Orthogonal to this flow we detect a second, 
curving flow, knots B 7-1 to B 7-3 . Neither flow seems to 
be associated with DR21(0H) itself: iLai. Girart fc CrutcheJ 



J2003I) and lVallee fc Fieed J2006r) both detect high- velocity 
CO emission associated with DR21(0H). However, the 
flow - mapped with a resolu tion of a few arcseconds by 
iLai. Girart fc Crutcherl ll2003l) - extends roughly east- west. 
Knot B 5-1 and possibly B 7-1 could be powered by this 
outflow. 

Jet B 7-1 to B 7-3 could instead be driven by the bright, 
nebulous near-IR source ~40" to the west of DR21(OH), la- 
beled DR21-IRS 6 in Fig|5| This source has Spitzer colours 
consistent with it being a reddened protostar, and is given 
reference letter "e" in Table |3] Sources "f" or "g" may 
likewise be associated with jet B 6-l:B 6-3, although the 
pre-main-sequence YSO ~10" east of B 6-2 is also well 
placed. Stars "f" and "g" have particularly steep spectral 
indices. Compari son of our images with the NH3 t hermal 
emission maps of lMangum. Wooten fc Mundvl lll992l) shows 



Observations of jets in DR21/W75 23 



r-i 













1 — 1 






M 











"3 


^ ~ ■ ■» 


1 










1 


































■ — ■ — m tmm m — ■ — ' 


1 










4 


tm 




.■" !!> 



IJ-E 



d v. 



*■ -B 1 



n 



4 



□ ■n 




r 




■ I 

"■ — p i ■ R ■ — ' 



Figure A2. Continuum-subtracted H2 l-OS(l) image of region "B" just to tlie north of DR21. The positions of YSOs are again marked 
with open circles and triangles; possible jet sources are given let ters. In th e uppermost contour plot, methanol masers associated with 
the DR21(OH) o utflow are marked with filled boxes IPlambeck fc Menten|[T99a: Kogan & Slvsh 1998), while water masers are marked 
with X symbols iForster et al ]|l978t iMangum. Wooten fc Mundvlll992^ . (See Fig^]for further details.) 



that these two jets he roughly midway between (and to the 
west of) two compact clusters of cores (labeled MM 1/MM 2 
and DR21(OH)N by this group). There is also no obvi- 
ous dust clump associated with sources "f", "g" or DR21- 
IRS 6/source "e", or with either jet, in the 850 /im map in 
Fig. El 

To the south in this figure we identify at least two 
(though possibly many more) coUimated jets radiating out- 
ward from the DR21 core and IR cluster: feature B 1-1 is 
probably a jet driving the knotty bow B 1-2:B 1-3, while 
knot B 4-1 is a possible bow heading northward from DR21- 
IRS 4 or a source in the vicinity of DR21 itself. The com- 



pact knot B 3-1 and the fainter features labeled B 3-2 are 
enveloped in faint, diffuse emission. These may delineate fur- 
ther coUimated flows extending northward from DR21. In- 
cluding the jets identified in the previous section, we find at 
least five coUimated jets radiating outward from the DR21 
850 /im peak and cluster of Hll regions and embedded IR 
sources (six outflows if we include the main DR21 outflow 
itseU). 

Continuing northward along the molecular ridge, we 
find a chain of bright H2 knots, B 9-1, which may be part 
of an east-west flow driven by the Spitzer YSO candidate 
source "h" at the western end of the B 9-1 chain. Alterna- 
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Figure A3. Continuum-subtracted H2 l-OS(l) image of region "C" centred on W75N (see Fig lAll for additional details). YSOs are again 
marked with open circles and trian gles; possible jet sources are assigned let ters. In the bottom-right con tour plot the positions of water 
masers are marked with X symbols iHunter et alJlOoS iTorrelles et al Il993); 2.7 mm continuum sources IShepherd. Testi fc Starkl2003l) 
are marked with pluses. Two of these mm sources, MM 1 and MM 5, coincide with W75N IRS1/W75N(B) and W75N IRS4/W75N(A), 
respectively. 



tively, this feature may represent a fragmented bow shock 
in a flow driven by an embedded source to the south. Pro- 
tostellar sources ERO 3, DR21-IRS 7 and DR21-IRS 8 are 
all possible candidates. 

There are also faint knots ^^20" south-east and ~20" 
north-east of DR21-IRS 10. Although this star is identified 
as a possible Class II or flat spectrum source in our Spitzer 
analysis (TableQl, it is also associated with a discrete source 
at 850 /xm (Fig.|Sj. Multiple unresolved sources may explain 
these seemingly contradictory data. IRS 10 certainly seems 
to be a strong source candidate for the coUimated flow as- 
sociated with B 11-1 and B 11-2. 

There are other faint H2 features distributed east of the 
ERO 3/FIR 1, FIR 2 and FIR 3 chain of dust cores, which 



we identify in the main panel in Fig. IA2I These are also evi- 
dent in Fig. 121 On morphological grounds alone, associating 
them with embedded sources along the molecular ridge is 
very difficult. Suffice to say that these features are generally 
within an arcminute (~0.9 pc) of the molecular ridge and 
may be part of east-west orientated jets. 



A3 W75N - region C 

In Fig. m we show a continuum-subtracted H2 image of 
W75N. The curving filaments of H2 emission associated with 
th e ends of the l arge-sc ale bipolar CO outflow, first identified 
bv lPavis et al.l 1^213) > S'l'^ o nce again evideii t (feat ures C 1- 
1 to C 1-4 and C 6-l:C 6-2). ITorrelles et alJ ^1997^ suggest 




Figure A4. Continuum-subtracted H2 l-OS(l) image of region "D" ~10' to the west-north-west of W75N. YSO positions are again 
marked with circles and triangles, although the Spitzer data are incomplete in this region; the dashed line marks the edge of the 3.6 fim 
and 5.8 /xm maps, while the dot-dashed line marks the edge of the 4.5 /^m and 8.0 /xm mosaics. (See Fig lAll for further details.) 



that this large-scale flow could be driven by VLA 1, based on 
the elongated structure of this source at radio wavelengths. 
The source is well aligned with the apices of the proposed 
bow shocks C 1-1 to C 1-4 and C 6-1 and C 6-2. Alterna- 
tively, knot C 1-5, within the wings of the C 1-1 to C 1- 
4 bow, may be part of a coUimated jet driving the bow 
shock, driven by the reddened K-band source W75N-IRS 2 
(Fig-EJ- Hydrogen recombination line emission has been de- 
tect ed from this star, so the source may be undergoing mass 
loss iMoore et alJll991lJ) . (The faint, diffuse emission a fur- 
ther arcminute to the south-west of knot C 1-2, noted by 
IShepherd. Testi fc Stark (2003), is also evident in the main 
panel in Fig. IA3I and as faint red features in Fig.|3) 

Knot C 7-1 to the north of the main W75N counter- 
flow lies along the axis of the extended knot C 4-1 and 
may be part of the same flow. A Spitzer YSO candidate 
and three dust cores traced in 2.7 mm continuum emission 



ijShepherd. Testi fc Starkl200,'jl all lie on approximately the 
same axis (although the millimetre peaks have no bright K- 
band nor 8.0 /im counterparts). There are regions of diffuse 
H2 emis sion between C 4-1 and C 7-1 (labeled NE-H and 
NE-I bv lShepherd. Testi fc Stari3 (l2003h ') which may also be 
part of this flow. Alternatively, this diffuse emission may be 
non-thermal, i.e. ambient gas fluorescently excited by neigh- 
bouring B stars. 

A further ~2' to the north, knots C 8-1 to C 8-5 trace a 
very well coUimated jet with almost the same position angle 
on the sky as the main W75N molecular outflow and the 
tentative C 4-1 to C 7-1 jet. However, if features C 9-1 and 
C 9-2 are also part of a coUimated jet, this flow is almost 
orthogonal to the other outflows. Two Spitzer YSO sources 
"j" and "k" are candidates for C 8-l:C 8-5, while source "1" 
likely powers C 9-l:C 9-2. In our K-band image "k" and "1" 
are both slightly extended, their conical nebulae opening 
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Figure A5. Continuum-subtracted H2 l-OS(l) image of L906E, region "E" (which is ~20' to the west of DR21). The positions of Spitzer 
YSO candidates are marked with circles and triangles; those that may drive the H2 jets are given letters. The dashed line marks the 
approximate edge of the four Spitzer mosaics. (See Fig lAll for additional details.) The white thumb-print-like patch that coincides with 
L906E-IRS sources 1, 2 and 3 is an artifact from the subtracted K-band image. It was produced by patches of dirt on the WFCAM field 
lens iDvc ct al 200&V 



in the directions of the C 8-l:C 8-5 and C 9-l:C 9-2 jets 
respectively. All three protostars have steep spectral indices 
(TablellJ, and together these three stars are associated with 
a compact molecular core, labeled W75N-N in the 850 /im 
map in Fig. |S|(see also Table |2j. 

We also mention knots C 5-1 and C 5-2 on the southern 
edge of the main W75N outflow. These features form part 
of a very faint, narrow structure that snakes back toward 
W75N-IRS 3. However, one could also argue that C 5-1 and 
C 2-2 form part of a coUimated jet, powered by source "i", 
which notably coincides with SCUBA core W75N-E (Ta- 
ble |21 although this peak may actually trace CO line emis- 
sion in SCUBA 850 nm filter bandpass, rather than dust 
continuum emission). 

Subtraction of the diffuse continuum that envelopes 
W75N(A) and (B) reveals a number of elongated H2 struc- 
tures. These are thought to be bona fide line-emission fea- 
tures, since they are also recovered in the H2 images of 



iDavis et al.l and ISheoherd. Testi fc StarM l)2003l) . 

The coUimated structure C 2-1 and knot C 2-2 may be part 
of another coUimated jet, although there is no bright K- 
band nor 8.0 fim source along this axis. Knots C 3-1 and 
C 3-3 may be bow shocks in the lobes of a bipolar jet, seen 
here as knotty feature C 3-2. Such a jet might be driven by 
a source midway between C 3-1 and C 3-3, although again 
no millimetre nor 8.0 fim peak is found here. 

Shepherd and collaborators devote a series of papers 
to examining the properties of the sources in the central 
W75N cluster, and to identifying the s ource or sources 
that drive the large-sca l e CO outflow dSheaherdl |200l|; 
[Shepherd. Testi fc Star3 l2003l : ISheoherd. Kurtz fc Testil 
!2004D . They identify a possible complex of multiple, over- 
lapping, small-scale outflows in W75N. They also point out 
that the VLA 1 jet lacks sufficient momentum to drive the 
large-scale CO outflow, and conclude that multiple sources 
may be producing the flow. Although our new near-IR 
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images support the idea of multiple flows in W75N (and its 
immediate surroundings), the data do not really shed any 
further light on this issue. 

A4 W75N west - region D 

Area D, to the west of W75N, was not observed in all four 
IRAC bands. Therefore, we can not use IRAC photometry 
to search for embedded YSOs. Even so, the region is of con- 
siderable interest. It includes a prominent bipolar molecular 
outflow labeled D 1-1 to D 1-8. On morphological grounds 
alone it seems likely that features D 1-1:D 1-2 and D 1-8 rep- 
resent bow shocks at the ends of bipolar flow lobes. Midway 
between these features there is a compact cluster of at least 
three near-IR sources. Two targets, W75N-IRS 10 and 11, 
are also 8.0 ^m point sources and are strong candidates for 
the outflow central engine (see also Fig. (KJ. W75N-IRS 11 
in particular has a [4.5]-[8.0] colour of ~3.1 (Table0 and is 
coincident with H2 emission. 

The elongated H2 "fingers" D 3-1 and D 3-2 in Fig. 1X1 
are reminiscent of the "bullets" in Orion (and perhaps knot 
B 4-1). These may originate from the bright K-band/8.0 fim 
sources W75N-IRS 7 (source "m") or W75N-IRS 8. The 
latter was only observed in two IRAC bands, though like 
W75N-IRS 7 it is quite red (Table 0. 

Further to the north, the H2 features D 4-1 to D 4-6 
again describe a coUimated jet. Features D 4-4 and D 4- 
5 are probably the wings of - and D 4-6 the cap of - a 
westward-moving bow shock. 

In Fig. 13 and |S] we label four other interesting K- 
band/S.O/xm sources, W75N-IRS 9, 12, 13 and 14. There 
is very faint II2 emission between sources W75N-IRS 9 and 
12, labeled D 11-1 in Fig. IA4I Additional faint H2 jets and 
knots are detected west of W75N-IRS 12 and around W75N- 
IRS 7/IRS 8. Sources 13 and 14 appear to be point sources 
within small ionised shells which we see clearly in Fig. 13 in 
scattered light and/or recombination-line emission. There 
is possibly some fluorescent H2 emission along the edges of 
each shell. We see no evidence for II2 jets associated with 
either source, however. 



cannot be used to search for the possible outflow source. 
However, they likely form part of a knotty jet. So might 
E 2-1 to E 2-3, which lie 20"-40" to the southeast of a pair 
of Class Il/flat spectrum sources. 

The faint knots E 6-1 and E 6-2 likely trace yet another 
coUimated flow, while feature E 7-1 appears to be a faint 
bow shock, possible driven by a source in the vicinity of 
L906E-IRS 3. 

Of the remaining H2 feature, E 1-1 is close to a nebu- 
lous K-band source L906E-IRS 6 (situated ~25' to the west- 
northwest of E 1-1). This source was just inside the 5.8 
and 8.0 /xm Spitzer field of view, and again appears to be 
fairly red (Table [TJ. 

This paper has been typeset from a T^niX/ I^TJ^X file prepared 
by the author. 



A5 L906E - region E 

Finally, to the west of DR21 there is a cluster of nebulous 
embedded sources and a bok-globule-like region of increased 
extinction in the near-IR. These features lie on the eastern 
edge of the Lynds dark cloud L906. Unfortunately this area 
is outside the field of view of the SC UBA 850 nm map, an d 
it was not discussed in great detail bv lMarston et al.H2004 . 
However, in Fig.|S|and Fig. lA5l one easily sees the rich nature 
of this region. 

We clearly identify a number of coUimated jets in area 
E. The wishbone-shaped feature E 4-1, E 4-2 and E 5-1 
to 5-6 is associated with a bright 8/im peak and spherical 
core of mid-IR emission, which we dub L906E-IRS 1. This 
source has a particularly steep spectral index and certainly 
qualifies as an ERO (Table 3J- Even so, the H2 knots are 
probably associated with independent flows, rather than a 
single wide-angled wind. 

The H2 features E 3-1 to E 3-4 are unfortunately beyond 
the edges of all four of the Spitzer IRAC images, so these 



